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This  dissertation  presents  the  results  of  several  types  of 
measurements  made  on  the  decametric  electromagnetic  radiation  from 
the  planet  Jupiter  using  a  very- long-baseline  interferometer.  Jovian 
signals  were  received  at  stations  located  at  Old  Town,  Florida, 
St.  Petersburg,  Florida,  Maipu,  Chile,  and  Bowling  Green,  Kentucky. 
The  interferometer  baselines  ranged  from  11,150  to  456,000  wavelengths 
of  the  18  MHz  receiving  frequency.  Interference  of  the  signals  was 
accomplished  by  using  the  cross  correlation  function  with  digitized 
data  in  a  program  for  the  IBM  360/65  at  the  University  of  Florida 
Computing  Center. 

By  interpreting  the  amplitude  of  the  cross  correlation  function 
as  the  fringe  visibility,  an  upper  limit  on  the  angular  size  of  the 
source  of  the  radiation  was  determined  to  be  0.1  seconds  of  arc,  if  it 
were  assumed  that  the  source  was  incoherent  and  had  a  Gaussian  brightness 
distribution.  The  high  correlation  of  the  envelopes  of  the  two  L-Bursts 
studied  using  the  very  long  baseline  (over  7,000  km)  indicated  that  the 
time  dependence  of  the  strength  of  the  burst  was  probably  intrinsic  to 
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the  source.  The  stability  of  the  fringes  formed  by  interference  of 
the  L-Burst  data  indicates  that  the  source  did  not  "jump  about"  more 
than  85  km  for  the  2.2  seconds' duration  of  the  1971  burst.  The 
"jumping  about"  motion  is  superimposed  upon  a  constant  drift  of  the 
source  with  respect  to  the  interferometer  power  pattern.  The  analysis 
of  a  series  of  S-Bursts  received  in  1970  indicates  that,  if  the  radia- 
tion from  Jupiter  is  beamed  and  the  beams  rotate  with  respect  to  the 
Earth- Jupiter  line,  then  the  minimum  sweep  rate  of  the  beams  must  be 
greater  than  10°/second  in  the  north-south  direction  and  greater  than 
0?4/second  in  the  east-west  direction.  The  method  of  sweep-rate 
measurement  involved  comparing  the  time  of  arrival  of  an  identifiable 
phase  point  on  the  heterodyned  Jupiter  signal  with  that  of  a  similarly 
identifiable  amplitude  point  on  the  envelope  of  the  detected  signal 
at  various  pairs  of  stations. 

While  the  use  of  the  interferometer  for  the  determination  of  source 
size  is  restricted  to  a  source  that  is  an  incoherent  radiator,  the 
remaining  analyses  place  new  limits  on  any  type  of  source. 
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CHAPTER  I 


A  GENERAL  INTRODUCTION 


1-1.  Jupiter's  Electromagnetic  Radiation  and  the  Wavelength 
Dependence  of  Its  Characteristics 

The  electromagnetic  radiation  received  from  Jupiter  is  the  main 
key  available  for  man  to  use  in  determining  the  nature  of  the  giant 
planet.  The  visible  spectrum  from  the  planet  is  that  of  sunlight 
reflected  off  the  tops  of  clouds.  The  presence  of  absorption  lines 
due  to  gases  above  the  reflecting  level  serves  to  modify  the  solar 
spectrum.   In  the  infrared  a  disk  temperature  of  about  130°K  was  measured 
by  D.  H.  Menzel  and  others  as  early  as  1926.  This  measurement  is  thought 
to  represent  the  temperature  at  the  cloud  tops,  or  even  somewhat  higher, 
if  the  radiation  is  considered  to  be  due  to  a  source  obeying  the  black 
body  radiation  laws.   It  was  reported  by  Mayer,  McCullough  and 
Sloanaker  (1958)  that  radiometer  measurements  at  a  wavelength  of  3  cm 
gave  a  disk  temperature  of  about  140°K.  The  value  is  in  approximate 
agreement  with  the  infrared  value,  and  it  is  thought  to  indicate  that 
the  radiation  at  this  wavelength  is  still  thermally  generated.  A 
measurement  at  10  cm  by  Sloanaker  (1959)  was  quite  different.  If  the 
source  of  the  radiation  received  at  this  wavelength  is  considered  to  be 
a  black  body,  then  its  temperature  must  have  been  about  600°K. 


Succeeding  measurements  at  the  longer  wavelengths  of  21  cm  and  31  cm 
seemed  to  confirm  the  trend  with  disk  temperatures  of  2000°K  and  5000°K, 
respectively.  Radhakrishnan  and  Roberts  (1961)  reported  that  the  source 
region  for  this  wavelength  was  about  three  times  the  diameter  of  the 
visible  planet  in  the  plane  of  the  equator  and  that  the  radiation  had 
a  moderate  amount  of  linear  polarization.  A  model  of  the  source  of  the 
radiation  at  this  wavelength  was  proposed.  The  model  accounts  for  the 
linear  polarization  by  having  relativistic  electrons  emit  synchrotron 
radiation  as  they  spiral  in  a  magnetic  field. 

While  radiation  between  the  frequencies  of  400  MHz  and  40  MHz  has 
yet  to  be  detected  with  certainty,  Jupiter  is  very  active  in  the  decameter 
wavelength  region.  The  decameter  region  contains  frequencies  from  40  MHz 
downward  to  about  5  MHz  where  the  earth's  ionosphere  becomes  opaque.  The 
flux  density  reported  by  Carr,  et^  al .  (1964)  increases  very  rapidly  with 
wavelength,  seemingly  with  a  spectral  index  on  the  order  of  +  8.  Radiation 
in  the  decameter  region  was  first  reported  by  Burke  and  Franklin  (1955) 
and  confirmed  by  Shain,  who  used  old  records  dating  back  to  1950.  Shain 
(1955  and  1956)  reported  that  the  source  of  the  radiation  did  not  emit 
isotropically  but  seemed  to  radiate  into  preferred  directions  in  a 
coordinate  system  attached  to  the  planet.  System  II  coordinates,  with 
a  rotation  period  of  9  55m40?6  (based  on  the  temperate  zone  visual 
markings),  provided  a  rotating  frame  in  which  the  radio  sources  appeared 
to  be  at  rest.  This  later  was  modified,  due  to  the  sources  drifting  in 
System  II,  to  a  longitude  coordinate  frame  with  a  rotation  period  of 
9h55m29s37  based  on  the  work  of  Shain  (1956),  Carr,  Smith,  Bollhagen,  Six  and 
Chatterton  (1961)  and  other  groups. 


1-2.  The  Nature  of  the  Decametric  Radiation 

The  most  noticeable  feature  of  Jupiter's  decametric  radiation  is  its 
extreme  intensity  variation  with  time.  Radiation  at  this  wavelength  comes 
in  the  form  of  noise  storms  that  last  on  the  order  of  a  few  minutes  to 
a  few  hours.   Gallet  (1961)  further  breaks  down  the  noise  storms  into 
burst  groups  (with  durations  on  the  order  of  minutes)  and  into  two  types 
of  pulses  that  make  up  the  burst  groups.  The  L-pulse  is  characterized 
by  having  a  length  of  from  0.1  to  5  seconds  while  the  S-pulse  is  shorter 
with  a  length  of  from  less  than  10  to  about  100  milliseconds.  Flagg  and 
Carr  (1967)  reported  that  bursts  substantially  shorter  than  this  have 
been  detected.  Theories  for  the  formation  of  the  two  types  of  pulses  are 
given  later  in  this  chapter. 

Due  to  the  extreme  time  variability  of  the  radiation,  it  is  usually 
more  convenient  to  use  the  probability  of  occurrence  of  an  event  as  the 
dependent  variable.  The  location  of  the  sources  with  respect  to  System 
III  was  accomplished  by  plotting  the  probability  of  occurrence  of  Jupiter 
radiation  versus  the  central  meridian  longitude  (that  longitude  in 
System  III  that  is  directly  below  the  earth  as  viewed  from  Jupiter).  Source 
A  is  located  at  approximately  250°,  Source  B  at  about  160°  and  Source  C 
at  about  320°  in  System  III.  The  source  regions  are  each  about  60°  wide. 
Bigg  (1964)  found  a  correlation  between  the  superior  geocentric  longitude 
of  Io,  the  innermost  Galilean  satellite,  and  the  probability  of 
occurrence  of  radiation  from  the  previously  described  sources.  Prediction 
of  Io-Source  events  can  now  be  made.  The  events  described  in  this 
dissertation  occurred  during  Io-B  predicted  storms.  The  Io-B  storms 
occur  when  Io  is  about  90°  from  superior  geocentric  conjunction  and 
Source  B  is  crossing  the  central  meridian. 


Jupiter  radiation  is  generally  of  relatively  wide  bandwidth.  Work 
by  Carr,  et  a[.  (1964),  Riihimaa  (1964)  and  Dulk  (1965)  characterizes 
the  bursts  as  having  bandwidths  on  the  order  of  a  MHz  with  center 
frequencies  that  move  in  time.  Some  interesting  features  in  the  form  of 
modulation  lanes  visible  in  the  amplitude- frequency  domain  have  been 
reported  by  Riihimaa  (1970).  Since  the  receivers  in  the  experiments 
described  here  have  bandwidths  on  the  order  of  a  few  kHz,  the  inter- 
action of  a  frequency-drifting  band  of  noise  with  the  narrow-band, 
fixed- frequency  receivers  must  be  borne  in  mind. 

Decametric  radiation  from  Jupiter  exhibits  strong  circular 
polarization.  Emission  from  Source  B  is  almost  always  polarized  in  the 
right  hand  sense  at  the  commonly  studied  frequencies  in  the  decametric 
part  of  the  spectrum. 

1-3.  A  Brief  History  of  Jovian  Decametric  Measurements 
Using  Long- Baseline- Interferometry. 

It  was  indicated  early  in  the  investigation  of  Jovian  radiation 
that  the  source  was  in  the  vicinity  of  the  planet  and  was  probably 
rather  small.  It  was  therefore  necessary  to  resort  to  interferometric 
methods  to  measure  the  size  of  the  source.  Gardner  and  Shain  (1958) 
found  no  burst-for-burst  correlation  between  signals  received  at  the 
ends  of  a  20  km  baseline.  Douglas  and  Smith  (1961)  reported  a  few  bursts, 
especially  of  the  S-type,  that  were  correlated  over  a  baseline  of  100  km. 
Slee  and  Higgins  (1963)  used  a  radio- linked  phc^e  sensitive  interferometer 
with  a  baseline  of  32  km  and  frequency  of  19.7  MHz.  They  found  that  the 
fringe  visibility  was  not  reduced,  hence  the  source  was  not  resolved, 
with  this  baseline. 


The  University  of  Florida  group  used  a  Brown  and  Twiss  Post- 
Detection-Correlation  Interferometer  before  1965  with  a  baseline  of 
55  km,  as  reported  by  Carr,  et  al_.  (1965).  May  (1965)  employed  data 
from  this  interferometer  to  determine  a  maximum  source  size  of  50  seconds 
of  arc.  A  phase  interferometer  was  also  used  by  this  group  that  employed 
the  cross  correlation  of  the  heterodyned  radio  frequency  noise.  Brown, 
et  al_.  (1969)  used  this  system  with  a  baseline  of  880  km  to  reduce  the 
maximum  possible  source  size  to  less  than  one  second  of  arc.  The  phase 
interferometer  of  the  University  of  Florida  group  used  tape  recorders 
for  the  storage  of  data  before  cross  correlation  was  done. 

Meanwhile  Dulk,  Rayhrer,  and  Lawrence  (1967)  used  detected,  tape- 
recorded  receiver  output  to  form  an  interferometer  of  the  Brown  and 
Twiss  type  with  a  baseline  of  175  km  at  34  MHz.  They  obtained  an  upper 
limit  of  3  seconds  of  arc  on  the  size  of  the  sources  of  both  L-  and 
S-bursts.  Dulk  (1970)  and  Stannard,  et  al.  (1970)  extended  the  baseline 
to  487,000  wavelengths  at  the  operating  frequency  of  34  MHz.  By  using 
the  instrument  in  the  phase  interferometer  mode,  they  were  able  to 
explore  the  phase  stability  of  the  source.  If  the  source  were  assumed 
to  be  incoherent,  they  found  that  its  maximum  size  was  less  than  0.1 
second  of  arc. 

A  preliminary  report  of  part  of  the  work  detailed  in  this  disser- 
tation was  made  by  Carr,  et  <d.  (1970)  in  the  form  of  a  history  of  the 
University  of  Florida  18  MHz  interferometry  experiment.  At  that  time  a 
456,000  wavelength  baseline  was  in  use,  but  only  small  samples  of  the 
data  had  been  hand  reduced. 


1-4.1.  Early  Theories  Concerning  the  Production  of  Jovian  Noise  Storms 

The  earliest  suggestions  as  to  the  origin  of  the  Jovian  decametric 
radiation  included  Jovian  lightning  storms  and  dynamo  effects  created 
by  slippage  of  the  cloud  bands  relative  to  one  another.  These  early 
proposals  ran  into  energetic  problems  in  accounting  for  the  high 
intensity  of  the  bursts. 

Warwick  (1963)  proposed  Cerenkov  radiation  of  relativistic 
electrons  as  a  source.  This  was  followed  by  Zhelesznyakov  (1965)  with 
a  theory  involving  plasma  waves  in  Jupiter's  ionosphere.  Ellis  (1965) 
is  the  first  to  have  published  a  model  that  involved  cyclotron  radiation 
from  electron  streams.  This  model  stands  as  the  direct  ancestor  to  most 
of  the  seriously  considered  theories  in  use  today.  The  discovery  of 
the  Io  effect  by  Bigg  in  1964  required  the  inclusion  of  an  interaction 
mechanism  between  the  satellite  and  the  emitting  region.  Gledhill 

(1967)  proposed  that  Io  interacted  with  an  implausibly  dense  plasma 
discus  about  the  planet  by  producing  waves  in  it.  Piddington  and  Drake 

(1968)  had  Io  disturbing  the  plasma  as  a  result  of  its  having  a  high 
degree  of  conductivity  or,  alternatively,  a  permanent  magnetic  moment. 
References  to  these  papers  are  contained  in  Schatten  and  Ness  (1971) . 

Goldreich  and  Lynden-Bell 's  theory  (1969),  which  is  discussed  in  the 
next  section,  requires  that  Io  be  about  100  times  as  conductive  as  the 
earth's  moon.  This  problem  is  alleviated  in  two  more  recent  theories. 
One,  by  Schmahl  (1970),  treats  Io  as  the  source  of  Alfvenic  waves  that 
travel  down  the  magnetic  field  lines  to  the  ionosphere  of  the  planet. 
The  waves  become  shockfronts  that  cause  radiation  from  coherent  electron 
groups.  The  theory  is  still  not  finished  in  that  it  does  not  develop 
the  interaction  mechanism  between  the  Alfvenic  waves  and  the  electromagnetic 
waves.  Another  more  recent  theory  by  Schatten  and  Ness  (1971)  does 


not  require  the  high  conductivity  for  Io  that  the  theory  of  Goldreich  and 
Lynden-Bell  does.  They  propose  that  Io  is  Moon-like  and  that  Jupiter's 
magnetic  field  lines  slide  through  Io  relatively  unaffected.  It  causes 
local  perturbations  in  the  pitch  angle  of  the  electrons  spiralling  in 
the  field  near  it,  however,  thus  causing  them  to  emit  radiation  when 
they  reach  the  top  of  the  Jovian  ionosphere.  Schatten  and  Ness  propose 
that  radiation  is  emitted  into  a  plane  perpendicular  to  the  flux  line 
that  passes  through  Io.  Their  paper  presents  a  computer  solution  to  the 
geometry  problem  involving  the  orientation  of  the  magnetic  dipole  of 
Jupiter  and  the  Earth-Jupiter  line.  The  solutions  agree  with  the 
locations  of  the  three  commonly  observed  sources,  and  also  the  rarely 
seen  fourth  source.  The  geometry  of  the  theory  shows  promise,  but;  as 
in  Schmahl's  theory,  the  mechanism  for  the  production  of  radiation  is 
not  fully  explored. 

1-4.2.  The  Theory  of  Goldreich  and  Lynden-Bell 

A  possible  explanation  of  the  effect  of  Io  on  the  Jovian  radiation 
is  given  by  P.  Goldreich  and  D.  Lynden-Bell  (1969).  They  propose  that 
the  satellite  acts  as  a  unipolar  inductor  which  generates  currents  along 
a  magnetic  flux  tube  connecting  it  with  Jupiter's  ionosphere.  The 
electron  currents  are  thought  to  undergo  cooperation  instabilities  near 
the  planet,  where  they  radiate  by  a  form  of  MASER  action  into  preferred 
directions  at  frequencies  just  above  the  local  gyrofrequency. 

From  the  mapping  of  the  source  of  the  decimeter  radiation  and  the 
investigation  of  its  linearly  polarized  component,  it  has  been  found  that 
a  magnetic  field  having  an  approximately  dipolar  geometry  can  be  assigned 


to  the  planet.  The  N  magnetic  pole  is  located  in  the  northern  hemisphere 
and  is  tilted  10°  from  the  planet's  rotation  axis  toward  the  System  III 
longitude  of  200°.  Figure  I . la  illustrates  the  geometrical  relations 
during  Io- related  storms. 

Goldreich  and  Lynden-Bell,  by  assuming  the  conductivity  of  Io  to 
be  about  that  of  silicates,  10"5mho/cm,  show  that  the  Jovian  magnetic 
field  would  diffuse  into  Io  in  about  one  day.  The  field  lines  passing 
through  Io  would  rotate  about  the  planet  at  Io's  orbital  speed.   If  Io 
has  a  finite  resistance,  the  flux  tube  will  not  be  entirely  trapped  but 
will  slowly  drift  through  it.  The  slippage  time  must  be  long  compared 
with  the  time  for  adjustments  to  be  made  in  the  rest  of  the  "circuit." 
The  adjustments  are  propagated  as  Alfven  waves.  Hence,  the  slippage  of 
a  flux  line  through  Io  must  take  longer  than  the  time  for  an  Alfven  wave 
to  travel  along  the  flux  tube  to  the  planet  and  return.  Since,  in  Io's 
reference  frame,  the  parallel  component  of  the  electric  field  must  vanish, 
an  electron  current  flows  toward  the  planet  along  the  outside  of  the  flux 
tube  and  returns  to  Io  along  the  inside  (the  side  nearest  Jupiter) .   In 
Jupiter's  frame,  there  is  an  electric  field  in  the  ionosphere  at  the 
feet  of  the  flux  tube  which  is  directed  toward  the  north  for  the  one  in 
the  northern  hemisphere.  The  feet  of  the  flux  tube  slip  relative  to  the 
rotation  of  the  planet  since  the  tube  rotates  about  Jupiter  at  Io's 
orbital  rate  (1.77  days).   In  the  ionosphere  the  flux  tube  has  an 
elliptical  cross  section.  When  the  ionization  of  the  ionosphere  is  high 
the  feet  of  the  flux  tube  will  be  pulled  ahead  of  the  satellite  by 
the  rotation  of  the  planet.  The  authors  of  the  theory  calculate  that  the 
feet  will  lead  Io  by  12°  under  these  conditions.  When  the  ionization  is 
low,  e.g.,  night  time  at  the  foot  of  the  flux  tube,  there  will  be  no 
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Figure  1.1  --  Geometrical  Relationship  of  Jupiter,  Io  and 
the  Earth  for  Io-Related  Storms. 
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dragging,  and  the  flux  tube  will  terminate  on  the  same  Jovian  longitude 
for  which  Io  is  on  the  Meridian.  During  the  transition  from  night  to 
day,  the  foot  of  the  flux  tube  will  be  swept  from  its  night  position  to 
a  point  some  12°  ahead  in  about  20  minutes  (i.e.,  at  Jupiter's  rotation 
rate  of  about  36°/hour) . 

The  weakly  relativistic  electrons  spiralling  in  the  flux  tube  can 
serve  as  negative  absorbers  if  their  momentum  distribution  is  sufficiently 
inverted.  The  greatest  amplification  occurs  in  the  extraordinary  mode 
for  propagation  at  right  angles  to  the  magnetic  field  as  seen  in  the 
coordinate  system  in  which  the  average  electron  streaming  velocity  is 
zero.  The  linear  Doppler  effect  causes  the  frequency  of  the  radiation 
to  be  raised  above  the  local  electron  gyrofrequency  as  seen  in  the  Jovian 
frame.   It  also  causes  the  radiation  to  be  directed  at  large  angles  with 
respect  to  the  electron  stream  velocity. 

The  radiation  from  10  keV  electrons  spiralling  up  from  the  planet 
along  the  flux  tube  will  be  directed  into  a  conical  surface  with  a 
half-angle  of  79°  that  opens  symmetrically  about  the  flux  tube  outward 
from  Jupiter.  The  apex  of  the  cone  is  located  at  the  foot  of  the  flux 
tube  in  Jupiter's  ionosphere.  Electrons  spiralling  toward  the  planet 
will  also  emit  radiation  into  a  similar  cone  which,  upon  reflection 
of  the  radiation  by  the  ionosphere,  will  open  outward  from  the  planet. 
Such  a  cone  will  intersect  Jupiter's  equatorial  plane  in  a  pair  of 
diverging  lines  whose  half-angle  is  between  65°  and  75°,  depending  on 
the  tilt  of  the  magnetic  axis  with  respect  to  the  Jupiter-Io  line.  The 
authors  of  the  theory  further  require  that  the  radiation  be  beamed  only 
into  selected  regions  of  the  conical  surface.  The  actual  degree  of 
beaming  anisotropy  cannot  be  predicted;  however,  for  reasons  which  seem 
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plausible,  the  assumption  is  made  that  the  preferred  beaming  angle  is 
at  13°  with  respect  to  the  east-West  direction  on  Jupiter.  The  assumed 
parameters  give  a  model  for  the  Source  B  emission  geometry  as  shown  in 
Figure  I. lb. 

Observations  show  that  the  sources  appear  to  be  narrower  at  higher 
monitoring  frequencies.   In  the  theory,  the  local  gyrofrequency  increases 
nearer  the  planet  such  that  the  source  of  the  higher  frequency  radiation 
must  be  closer  to  Jupiter.  The  axis  of  the  radiation  cone  is  tilted 
farther  from  the  equatorial  plane  since  the  flux  tube  is  more  nearly 
parallel  to  the  magnetic  dipole  at  lower  altitudes.  The  greater  the 
tilt  of  the  cone  axis  the  smaller  the  angle  is  between  the  lines  of  inter- 
section between  the  cone  and  the  equatorial  plane.  This  makes  the  source 
seem  narrower  when  viewed  from  Earth  since  radiation  reaching  the  earth 
can  travel  only  along  the  lines  of  intersection. 

The  theory  gives  a  model  for  the  formation  of  the  S-Burst.  It 
proposes  that  a  column  of  electrons  spiralling  on  the  surface  of  the  flux 
tube  emits  coherently  for  a  period  of  time.  The  time  duration  of  the 
S-Burst  is  determined  by  the  coherence  time  of  the  column  of  electrons. 
The  bandwidth  of  the  radiation  is  determined  by  the  altitude  range,  hence 
the  range  of  local  gyrofrequencies ,  spanned  by  the  column.  S-Bursts  may 
be  formed  by  electrons  moving  in  either  direction  with  respect  to  the 
planet  during  the  time  for  emission  of  radiation. 

The  source  of  the  radiation  energy  is  the  transverse  kinetic  energy 
of  the  keV  electrons  spiralling  in  the  magnetic  field.  The  electrons  have 
undergone  acceleration  from  an  energy  of  about  4  eV  to  the  keV  region 
by  the  electric  field  generated  by  Io  moving  in  Jupiter's  magnetic  field. 
The  ultimate  source  of  energy  for  the  radiation  resides  in  the  orbital 
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energy  of  Io.  Using  plausible  values  for  the  number  density  of 
radiating  electrons,  etc.,  the  authors  show  that  the  doubling  of  Io's 
angular  momentum  is  comfortably  long  compared  with  the  lifetime  of  the 
solar  system. 

A  serious  weakness  in  the  theory  occurs  in  the  assumed  value  for 
the  conductivity  of  Io.  It  was  necessary  to  assume  a  conductivity  of 
about  two  orders  of  magnitude  greater  than  that  of  the  moon,  an  assumption 
which  may  not  be  realistic. 

We  wish  to  consider  the  theory  of  Goldreich  and  Lynden-Bell  in  this 
dissertation  since  it  predicts  the  beaming  of  the  decametric  radiation. 
If  the  beams  sweep  across  Earth,  they  should  arrive  at  one  station  of  an 
interferometer,  such  as  ours,  before  they  arrive  at  another  station.  We 
will  report  on  the  measurements  of  arrival  times  in  a  later  chapter. 

1-4.3.  The  Theory  of  Douglas  and  Smith  for  the  Formation  of  L-Bursts 
by  the  Solar  Wind 

J.  N.  Douglas  and  H.  J.  Smith  (1967)  proposed  a  theory  to  account 
for  the  L-burst  structure  in  the  Jovian  radiation.  Two  types  of 
experiments,  both  involving  spaced  receivers,  led  to  their  conclusion. 
Gardner  and  Shain  (1958)  had  found  that  the  burst  envelope  correlation 
was  often,  but  not  always,  poor  over  baselines  of  only  a  few  tens  of 
kilometers.  Carr,  et  al .  (1964)  obtained  a  similar  result  over  a  base- 
line of  7000  km.  The  cases  of  positive  correlation  indicated  that  the 
cause  of  the  amplitude  variation  in  the  burst  must  be  more  distant  than 
the  Earth's  ionosphere,  since  such  local  variations  should  not  influence 
the  signals  at  both  receivers.  Douglas  and  Smith  (1961)  started  a 
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continuing  program  of  monitoring  Jovian  radiation  for  the  purpose  of 
trying  to  detect  differences  in  arrival  times  of  a  given  burst  at 
several  stations.  Results  of  this  work  (Douglas  and  Smith,  1967) 
indicated  that  there  is  a  marked  dependence  of  the  delay  time  and  the 
order  of  arrival  at  receivers  in  the  east-west  orientation  on  the 
number  of  days  before  (or  after)  opposition  of  Jupiter. 

In  their  paper  (1967) ,  they  argue  for  a  random  distribution  of 
phase-changing  irregularities  in  the  solar  wind  located  in  the  space 
between  Jupiter  and  the  earth.  They  further  require  that  these 
irregularities  be  remote  from  both  the  earth  and  Jupiter.  The 
irregularities  are  stable  in  structure.  Their  motion  causes  an 
isophotal  pattern  to  drift  past  the  earth-based  receivers  at  a  rate  of 
several  hundred  kilometers  per  second.  Mariner  II  data  concerning  the 
solar  wind  indicated  that  the  velocity  of  the  wind  did  not  change 
appreciably  with  distance  from  the  sun.  Douglas  and  Smith  proposed  that 
an  inhomogeneity  in  the  solar  wind  would  cause  phase  changes  in  the 
radiation,  creating  a  drifting  isophotal  pattern  as  seen  at  the  fearth. 
According  to  their  calculations,  using  data  extrapolated  to  the  region 
of  the  Earth-Jupiter  line-of-sight,  they  were  able  to  account  for  the 
drift  rate  of  the  isophotal  pattern.   Furthermore,  their  theory 
predicted  that  before  Jovian  opposition  the  isophotal  pattern  should 
drift  eastward,  and  after  opposition  the  pattern  would  drift  westward. 
The  drift  direction  is  determined  by  the  sense  of  the  solar  wind  velocity 
vector  projected  on  a  plane  perpendicular  to  the  line-of-sight  from  the 
earth  to  Jupiter.   It  is  clear  that  the  vector  will  switch  directions 
at  opposition  if  the  solar  wind  is  assumed  to  move  in  a  radial  direction 
from  the  sun.  The  predicted  drift  direction  was  shown  to  agree  with 
experimental  data. 
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1-5.  The  Purpose  of  the  Experiments  Described  in 
This  Dissertation 

Theories  are  the  intellectual  creations  of  mankind  in  the  same 
sense  as  are  the  more  conventional  arts.  In  the  natural  sciences, 
however,  man  has  added  the  constraint  that  his  theories  must  represent 
the  way  in  which  a  class  of  natural  phenomena  may  be  caused.  There  is 
thus  the  requirement  that  any  prediction  concerning  a  natural  phenomenon 
made  by  a  theory  must  agree  with  the  measurements  made  on  that  natural 
phenomenon.   In  case  of  disagreement,  difficulties  may  lie  both  in  the 
theory  and  in  the  interpretation  of  the  measurements.  It  is  generally 
not  possible  to  make  THE  measurement  that  will  unequivocably  make  the 
choice  between  two  or  more  theories  obvious. 

In  astronomy,  man  is  faced  with  taking  data  from  phenomena 
that  are  secondarily  related  to  the  information  he  is  seeking.  It  is 
therefore  necessary  to  examine  not  only  the  performance  of  the  theory  in 
question,  but  also  to  be  aware  of  possible  failures  of  theories  concerning 
the  way  the  information  is  transmitted.  A  case  in  point  is  the  deter- 
mination of  the  features  of  Jupiter's  magnetic  field.  The  radio 
astronomer  received  electromagnetic  radiation  that,  according  to 
Maxwell's  equations  and  the  theory  of  electrical  measurements,  could  be 
shown  to  be  partially  linearly  polarized.  The  synchrotron  process  had 
been  shown  to  be  a  source  of  linearly  polarized  radiation  having  the 
spectral  characteristics  observed  by  the  astronomer.  Since  this  demon- 
stration had  been  made  on  Earth,  one  had  to  make  the  assumption  that  such 
a  process  could  be  carried  on  elsewhere  given  similar  conditions.  This 
assumption  is  part  of  a  set  summarized  by  Newton  in  the  Principia  as  a 
list  of  Rules  of  Reasoning.  One  condition  for  the  occurrence  of  the 
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synchrotron  process  was  the  presence  of  a  magnetic  field.  Jupiter,  then, 
must  have  a  magnetic  field.  Further  measurements  of  the  distribution 
of  the  source  of  radiation  led  to  the  current  understanding  of  the 
geometry  of  the  field. 

At  the  same  time  that  an  inquiry  is  being  made  into  the  origin  of 
one  natural  phenomenon,  i.e.,  the  radiation  from  Jupiter  in  the 
previous  example,  many  other  phenomena  are  being  investigated.  There 
is,  for  example,  exploration  of  many  other  features  of  Jupiter  taking 
place,  such  as:  experimental  studies  of  the  atmospheric  features  of  the 
planet,  of  the  red  spot  and  theoretical  studies  of  the  interior  of  the 
planet.  The  aim  of  man's  inquiry  seems,  in  this  case,  to  be  that  he 
is  hoping  to  have  at  some  future  time  a  complete  picture  of  Jupiter 
which  can  be  added  to  the  rest  of  the  pictures  he  has  made  of  the 
universe.  By  that  time,  if  he  is  successful,  he  will  feel  at  home  in  a 
familiar  world,  a  world  of  his  own  intellectual  creation. 

The  experiments  in  this  dissertation  are  proposed  to  give  man  a 
little  more  insight  into  the  natural  phenomena  that  he  is  hoping  to 
understand. 

List  of  Experiments 

1.  The  angular  size  of  the  source  of  Jupiter's  decametric  radiation 
will  be  determined  using  both  a  phase  and  an  intensity  interferometer. 

2.  The  stability  of  the  position  of  the  source  will  be  examined 
using  the  fringe  pattern  generated  by  the  interferometer. 

3.  The  arrival  times  of  S-  and  L-Bursts  at  several  stations  will  be 
examined  to  provide  experimental  information  in  connection  with  the  theories 
of  Goldreich  and  Lynden-Bell  and  of  Douglas  and  Smith.  A  new  method  for 
eliminating  the  residual  timing  errors  will  be  used. 


CHAPTER  II 

GENERAL  INTERFEROMETER  THEORY 

The  theoretical  structure  for  the  understanding  of  the  phase  and 
intensity  interferometer  is  developed  in  the  first  two  sections  of  this 
chapter.  Since  the  interference  pattern  of  the  received  signals  from 
Jupiter  is  actually  obtained  from  the  computer,  the  last  section  of 
this  chapter  is  devoted  to  establishing  the  connection  between  the 
classical  forms  of  the  interferometer  and  the  experimental  technique 
used  here. 

II-l.  The  Phase  Interferometer 

Classically,  the  phase  interferometer  adds  the  signal  voltages  from 
two  spaced  antennas  at  some  central  point.  The  fluctuation  in  the 
combined  signal  power  is  dependent  on  the  size  and  location  of  the  source 
with  respect  to  the  interferometer  power  pattern.  Consider  an  inter- 
ferometer formed  by  placing  mixers  between  each  antenna  and  the  signal 
adding  point.  A  system  of  this  nature  allows  the  signals  to  be  stored 
on  magnetic  tape  and  the  interference  to  be  accomplished  under  more 
convenient  conditions . 

Figure  II. la  shows  a  block  diagram  of  the  system  to  be  considered. 
The  left  half  of  the  system  is  called  receiver  one  and  has  an  input 
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Diagram  of  the  Classical  Interferometer, 
Graph  of  the  Output  S  Versus  Beam  Angle  $  . 
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frequency  u>.,  a  local  oscillator  frequency  u>Lo,  and  a  low-passed  output 
frequency  a)Q.  The  right  side,  receiver  two,  has  similar  variables  that 
are  given  by  ul,  u^,  «J.  For  a  single  frequency  the  input  voltage  to 
receiver  one  is  E_  =  E^  (?)  cos  fa.t  +  <F),  where  ¥  is  the  phase  lag 

introduced  by  the  extra  path  length  y2  with  respect  to  the  point  halfway 
between  the  antennas.  Its  value  is  4-  =  ff  Cs+a)  sin  g  >  Thc  input  voltage 

A 

to  receiver  two  is  EJ  =  £•„(«  cos  fojt  -  <F),  where  f  has  the  same  value. 

The  amplitude  factors,  E^  and  E|q,  are  angle  dependent  since  each  antenna 

has  a  power  pattern  appropriate  to  its  design.  The  angle  variable,  ?  =  <j>0-<(., 

denotes  a  direction  in  the  far  field  pattern  of  the  interferometer 

measured  from  the  center  of  the  main  lobe.  The  local  oscillator  injection 

voltages  take  the  form:  E   =  E    cos  (a.  t)  and  E'  =  E'   cos  (o)«  t  +  n) , 

LO    L00      Lo       Lo    Lo0      Lo 

where  n  is  the  phase  at  time  =  zero.  Before  low-pass  filtering  the  outputs 

of  the  mixers  are  E  =  E   (?)  cos  (w.t  +  y)  Er   cos  (to  t)  and 

u    1U  *       Lo       Lo 

E0  =  Eio  (0  C0S  ^i*  ~  ^  ELo  C0S  ^lo*  +  n)>  Each  of  these  voltages   • 

contains  terms  that  oscillate  with  frequencies  to.  ±  ui   and  oj1  ±  u1 

*    Lo      i    Lo 
The  summing  point  will  receive  only  the  difference  frequencies  due  to  the 

low-pass  filters  in  each  mixer. 

Some  restrictions  are  now  placed  on  the  variables  that  are 

multiplied  in  the  mixers.  Let  oj.  =  u>.' .  If  the  source  emits  over  a  wide 

frequency  range,  these  equations  apply  to  each  frequency  in  the  bandwidth 

of  the  low-pass  filters.  Specifying  u*  =  a       +  y   provides  a  means  of 

Lo    Lo 

distinguishing  the  fringes  from  scintillations  in  the  interplanetary 
medium.  It  has  been  shown,  Dulk  (1970),  that  these  scintillations  have 
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a  power  spectrum  that  contains  very  little  power  above  a  frequency  of  a 
few  Hz.  If  y  is  greater  than  this  amount,  then  several  fringes  will 
pass  during  one  phase  fluctuation  due  to  the  interplanetary  medium.  It 
is  convenient  to  let  ELq  =  E'Lq  =  E  .  and, since  the  individual  antennas 


o      o 


are  steered,  E  (0)  =  E'  (0)  =  1. 

The  voltage  at  the  summing  point  is  given  by  E   =  E  +  E'  = 

E^cos  (oij.t  +  V)   cos  0LQt)  +  cos  (<uL  t  +  yt  +  n)  cos  (u.t  -  V)}.  Upon 

using  a  trigonometric  identity  and  the  condition  that  |u>  -  a     I  »  y 

i    Lo' 
the  summing  point  voltage  becomes  E   =  E  {cos  rfoi.  -  to  It  -  HI 

cos  (*  +  II  +  i  )}.  The  power  pattern,  after  smoothing  on  a  2^/(10.  -  u  ) 
^   *  i    Lo 

time  scale  and  normalizing,  is  P  (5,t,n)  =  cos2  (Y  +  II  +  H)  .  The  power 

n  2   2 

pattern  can  be  written  as  Pn(?,t)   =  1  +  cos  (6  +  yt) ,  where  6  =  2*  is 
the  total  phase  difference  due  to  path  length  caused  by  the  separation  of 
the  antennas.  The  y  term  gives  temporal  fringes  if  e  is  held  constant. 
The  phase,  n,  does  not  depend  on  the  source  and  can  be  chosen  to  be  zero 
for  convenience. 

The  observed  flux  density  from  a  source  is  given  by  S.(<|>  ,t)  = 

/  B(<J>)  pn(<l>0   -   <J0  d<j>,  where  B(<f>)  is  the  true  source  brightness  distri- 
source 

bution,  <|>0  is  the  hour  angle  of  the  main  lobe  of  the  interferometer  power 

pattern,  and  $   is  the  hour  angle  of  a  point  in  the  source.  Using  the 

previously  calculated  normalized  power  pattern,  P  (£,t),   the  observed 

a  a  " 

flux  density  becomes  S0(*0,t)  =  J2  B($)  d*  +  J2   B(<(.)  cos  (9  +  yt)  d<|)  where 

_«  a 

2  "7 

the  angular  width  of  the  source  is  a  and  5  =  <J>  -  <J>.  The  integral 


0 


s0  =  J  BC<}')  d<()  is  the  total  flux  density  from  the  source.  If  a  is 
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small  compared  with  the  lobe  spacing,  the  values  of  «J>  for  which  the 
second  integral  has  any  appreciable  size  are  limited  to  the  region 
I*  "  *„|  «  t.  Under  this  condition  9  becomes  6  =  2ir  s  (*„-$). 

a_ 

The  second  integral  takes  the  form  /2  B($)  cos  [2tt  (s  i>     +  vt)] 

o_  X      0 

"2 

a 


'  iH<PJ   sin   12tt    [s. 
a  A     o  a 


cos    (2ir  sx  <j>)   d<f>  +  /^  Bft)   sin   [2ir    (s^   *fl  +  vt)]   sin   {2v  S,    *)   d*,  where 

~2 


y   =  2™  and  sA  =  2S.  If  SQ  »0,»x.t)  is  written  as  SQ  (^V^  = 
s0  [1  +  V  (<j>0,sx,t)],  the  visibility  function  can  be  identified  as 


a 


V(*0,sA,t)   =    (i_)    {cos    [2tt    Cs^  <t>Q  +  vt)]  J2  Bfo)    cos    (2ir  s     $)  d*  + 

"2 

a_ 

sin   [2tt    (s     <J.0  +  vt)]   /2  BO)   sin    (2tt  S,    <j.)   d*}.      Let  VWn,s„t)  be 

ex  A  0   A 

"2 

written  as  VC^s^t)  =  VQ  (sx)  cos  [2ir  (sx  *q  -  sx  A  *Q)  +  Yt].   1) 
Figure  II. lb  is  a  graph  of  S  (+0,s  ,t)  plotted  against  $„  in  which  the 
local  oscillator  offset,  y,   is  made  to  be  zero. 

It  will  be  useful  for  later  derivations  to  identify  the  following 
terms  in  an  expansion  of  V  (<j>0,s  ,t)  above  (eq.  l) : 

A 

a 


V0  CV  C°S  C2U  SX  A  V  =  ~  /   BC<f))  COS  O  s   ♦)  d* 

S_   a  A 

0  'J 


*   (i")  F    BC*) 


i.g-/  F    BC«  2) 

°o   cos 


and 
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V  Cs  )   sin  (2ir  s  A  a  )  =  (i-)  /2  B(<J.)  sin  (2tt  s,  <J0  d<J) 
A  0   « 

'2 

+  %'  "sin  B»>  '  3> 

where  Fcog  and  F    are  the  Fourier  cosine  and  sine  transforms  of  the 
brightness  distribution  B(<J>).   In  order  to  make  the  identification  of 
the  transforms  it  is  necessary  to  require  that  there  are  no  other 
sources  of  appreciable  strength  in  the  antenna  pattern. 

If  Equation  3)  is  multiplied  by  i  and  added  to  Equation  2) ,  and 
the  Fourier  transform  is  taken,  the  following  equation  results: 

B(V  =  S0  C   V0  CS^  e"i27rWAV  dsx  .       If  the  visibility 

function's  amplitude,  V  (sx) ,  can  be  determined  for  all  baseline  lengths, 
the  brightness  distribution  of  the  source  can  be  reconstructed  by  an 
integration  over  the  baseline.  In  the  classical  interferometer  this 

S      S 
quantity  is  determined  by  evaluating  the  equation  V  (s  )  =  max    min 

0   *    S    +  S 

max    min 

for  many  different  baselines.  A  practical  limitation  on  this  method's 
use  is  that  many  measurements  are  required  with  increasingly  longer 
baselines.   In  the  case  of  a  small  source  emitting  at  long  wavelengths, 
as  for  Jupiter,  the  earth's  diameter  may  form  the  limit  on  the  integral. 
An  added  requirement  was  made  when  the  Fourier  transform  took  place.  The 
source  was  required  to  be  the  only  one  in  the  antenna  pattern.  This  is 
reasonably  true  of  Jupiter.  Subject  to  these  conditions,  the  brightness 
distribution  of  the  source  can  be  determined  from  a  knowledge  of  the 
visibility  function,  V  (<{>q,s.  ,t) . 


22 


It  is  possible  to  calculate  V  (<f>  ,s,  ,t)  by  a  method  that  is  more 

OX 

easily  realized  on  a  digital  computer.  The  technique  uses  the  cross 
correlation  function  which  is  defined  as 


E  (t)  E-  (t  +  t) 
C  (s,,C,t)  = 


[E02  CO]*  [E'2  (t  +  t)]* 

where  EQ  (t)  is  the  voltage  output,  after  low-pass  filtering,  of  the 
first  mixer  and  E^  (t  +  t)  is  that  of  the  second  mixer  at  time,  t,  later 
(see  Figure  II. 2a).  The  voltages  are  given  by 

a 

J 

E0    (t)   =  a  J     E  ar   (<j>)   cos   [{to     (<J0   -  w     }  t  +  n     (♦)  +  ¥     (♦)]   d<j>  and 
ot  —  *■  JjO  r  r 


E'    (t  +  x)  =  b  J     E  a.    (V)   cos   [{co.    (♦•)   -  to.    }(t  +  t) 
a  £  *  ko 

~2 

where  a  (<J>')  and  a  0)  are  the  amplitudes  of  the  frequencies  co  and  10  . 
%  r  £     r 

The  same  restrictions  as  for  the  classical  interferometer  are  added,  i.e., 
E.   =  E*   and  to'  =  to   +  y.  The  phases  of  each  component  frequency 

of  the  source  are  n  0)  and  n  O1),  and  the  initial  relative  phase  of 

the  local  oscillators  is  ignored.  It  should  be  noticed  that  the  amplitudes 

and  phases  of  the  £  and  r  components  depend  on  positions  in  a  random  way. 

It  is  further  required  that  B  (<J>)  a  E  a„  (cb)  a  U)  .  The  summation  on 

£,r  r     * 

H   and  r  is  carried  only  over  the  frequencies  passed  by  the  low-pass 
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filters  to  the  multiplier.  The  filters  are  assumed  to  have  a  flat 
response  inside  their  passband  with  no  response  outside,  and  the  band- 
width of  each  is  much  less  than  the  center  frequency. 

The  numerator  of  C  (s,,£,t)  becomes  E  ft)  E  ' ft  +  t)  = 

A  0  0  ' 


ab 


f"/t     L     J  ar     W   cos   l>ro   C*)   t  +  rir   (♦)   +  V     (*)]  d<j> 

/#l  J  a£  (<{.')   cos   [«ao    (<(.'  +  Y)(t  +  T)  +  n£   (+•)   -  ¥     (<(.')]  d<(.»  dt 
where  u>       O)   =  w     (*]   -  u       (♦)   and  »       (♦»)   =  *a    [*']    -  m.      (*»).     The 

x  LO  Jo  0  x,  LO 

variable  ¥  is  defined  as  before,  !  ■  i  s  ft,  -  ♦),  where  the  small 
angle  approximation  has  been  made.  The  integration  period  T  is  taken 
as  being  long  with  respect  to  2-n/u       and  short  compared  to  2tt/y  and  the 
fringe  rate  due  to  the  source  moving  across  the  interferometer  pattern. 
A  diagram  of  the  interferometer  is  given  in  Figure  II. 2a.  Such  a  system 
is  called  a  multiplying  interferometer  and  has  the  varying  part  of  the 
observed  flux  density  as  its  output.  By  rearrangement  of  the  integrations 


the  voltage  product  becomes  E  (t)  E  •  (t  +  t)  =  /  d<j>  —  f  f    z  a  f41 

<P       t  <p'  r,£ 

cos  [  u   (*)  t  +  n  (♦)  +  ¥  (*)]  ao  (♦')  cos  [(w   +  Y)(t  +  t)  + 

n£  C* ' 3  -  ¥£  (*')]  d<J>'dt.  It  is  necessary  to  examine  the  behavior  of  the 
amplitudes  and  phase  factors  under  the  integrations  over  time  and  <j>.  The 
amplitudes  and  phases,  a  O)  and  n  (*),  of  the  frequencies  are  random 
functions  of  position.  This  requires  that  the  integration  over  $'  be  zero 
except  where  $'  =  <J>.  Also  if  the  integration  over  time  at  any  point  is 
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Figure  II. 2  —  (a)     Diagram  of  the  Multiplying  Interferometer, 
(b)  ;    Graph  of  the  Output  Sfl  as  It  Depends  on  * 
(solid  line)  and  on  <frfl  (dashed  line). 
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considered,  the  integral  will  be  zero  unless  u>   -co   <  — .  This 

ro    £o   T 

is  caused  by  the  random  relation  of  phases  for  two  different  frequencies 

and  the  integration  time.  If  this  term  serves  as  the  integrand  for 

an  integral  over  position,  the  result  will  be  zero  unless  u   =  tu  . 

ro    £0 

which  is  again  due  to  the  random  nature  of  the  phases  and  amplitudes 
of  different  frequencies  with  position. 

We  are  left  with  the  product  term  having  the  form  E.  (t)  EQ ' (t  +  x)  = 

cib  ( 

Y~  J     z  ar  ft)  a.  (<J>)  cos  [£t  +  2¥  (<j>)]  d<j>,  where  the  integration  time 
r=£ 

is  short  compared  to  2tt/Y  and  to  the  time  for  the  source  to  pass  through 

one  fringe.  This  is  the  same  form  as  the  varying  component  of  the 

observed  flux  density  as  derived  for  the  voltage  adding  interferometer. 

The  denominator  of  C  (s  ,£,t)  can  be  written  as  [E2  (.t)]%   = 

f—  /  [  /  /  E  a*  (40  a  (*•)  cos  [a>  t  +  n  (♦)+¥]  cos 
T  r,s  ro     r       r 

[«sot  +  ng  (♦')  +  V  (♦•)]  d$»  d$]  dt}*5  ,  where  the  product  of  the  spatial 
integrals  has  been  replaced  by  the  double  integral  over  the  position 
variables.  By  the  arguments  used  in  the  previous  part,  $ ■  =  $  and 


"ro  =  u£0  or  the  inteSral  will  be  zero.  Then  [E  2  (t)]^  = 

2 

fj  d+  f-/   E  a  (+)  ae  (<{.)  cos2  [u  t  +  n  (*)  +  V  (<(.)]  dtK2.  After 
r=s                  ro     r       r 


the  time  integration  this  becomes  [E  2  (t)]2  =  {  —  /  Z  a  (<t>)  a  ((())  d*}"5, 

r=s 

where,  upon  substituting  B  (<J>)  for  the  sum  on  frequency,  we  recognize 
[EQ2  ft)]3*  =  {SLj   B  (*)  d^K2"  to  be  a  S^/ZI.  A  similar  term  for  the  rms 
value  of  the  second  receiver's  voltage  is  also  demonstrable. 
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The  cross  correlation  function  can  then  be  written  as 


re             ,               E0    (t)   E0      (t  +  T)                   f  V(sA,*0,t)S0 
C(sw?,t)  =    = i— ° 2. 


'X- 


[E  2    {t)]h   [E,2    (t  +  tjj%  a       ^  b       ^ 

°  TS„7S„ 

•2     0    V2     ° 


S0   Vo    (SX}    C0S    [Yt   +   27r  SA    ^0   "   A<V] 


s 

0 


It  is  then  seen  that  the  cross  correlation  function  is  just  the  fringe 
visibility.  The  amplitude  of  the  fringes  generated  by  plotting  the 
cross  correlation  function  against  time  is  just  V.  (s,).  The  cross 

0    A 

correlation  function  lends  itself  to  being  calculated  with  digitized 
data  in  the  computer.  Figure  II. 2b  is  a  graph  of  the  output  of  this 
form  of  interferometer  plotted  against  time.  The  envelope  represents 
fringes  due  to  the  motion  of  the  source  with  respect  to  the  interfero- 
meter fringes  and  the  solid  curve  is  the  cross  correlation  function. 
The  amplitude  variations  are  due  to  the  relative  local  oscillator  offset 
and  source  motion. 

II-2.  The  Intensity  Interferometer 

A  new  form  of  interferometer  was  proposed  by  R.  Hanbury  Brown  and 
R.  Q.  Twiss  (1954)  that  had  the  advantages  of  being  operable  over  a 
very  long  baseline  and  of  being  relatively  insensitive  to  phase  distortions 
caused  by  electron  density  variations  in  the  ionosphere.  A  diagram  of 
the  system  is  given  in  Figure  II. 3a.  The  Michelson  (phase)  interfero- 
meter is  sensitive  to  phase  shifts  that  occur  to  the  signals  before  they 
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Figure  II. 3  --  (a)      Diagram  of  the  Brown  and  Twiss  Post- 
Detection-Correlation  Interferometer,  (b) 
a  Sample  Record. 


/.a 


reach  the  summing  point.  The  system  described  here  is  not  as  sensitive 
to  this  problem  since  the  data  sent  from  the  remote  station  is  square- 
law  detected  and  filtered.  This  data  undergoes  two  types  of  treatment 
at  the  time  that  interference  is  applied.  Each  signal  is  linearly 
detected,  filtered  and  plotted  on  a  chart  recorder.  A  third  channel  is 
formed  by  multiplying  the  original  signals,  filtering  and  then  plotting 
the  product.  A  delay  network  is  incorporated  into  one  leg  to  compensate 
for  the  time  delay,  t,  caused  by  the  data  transmission  link.  The 
normalized  correlator  output  c  (£,t)  is  given  by  the  equation, 


\   Ct)  •  s2  (t  -  t0) 

c      =  <rvM]*-y<ivct-vi*-v  '  "here  the 


numerator  is  the  filtered  multiplier  output  signal  at  a  given  time,  t, 

t         1 
and  [Sj  ]2  and  [S^]"1  are  the  values  of  the  individual,  linearly  detected 

channels  at  the  same  time.  P   and  P   are  the  average  values  of  the 

Nl      N2 

low-frequency  filtered  outputs  observed  when  the  source  is  outside  the 
antenna  beam.  C  (£)  is  formed  by  averaging  c  (£,t)  over  a  number  of 
passes  of  the  source  through  the  beam.  The  variable  S,   is  the  separation 
of  the  two  antennas. 

Brown  and  Twiss  show,  in  the  first  part  of  their  paper,  that  a  term 
they  call  the  correlation  coefficient, p  (£) ,  is  given  by 

[F   2  (s.)  +  F  .  2  (s  )]% 

L  cos  v  A'    sin  v  A  J    ,. 

p  (£)  =      ~~ for  the  Michelson  interferometer. 

F    CO) 
cos 

It  is  seen  that  the  numerator  can  be  formed  by  referring  to  the  square 
root  of  the  sum  of  the  squares  of  Equations  2)  and  3)  in  the  previous 
section.  The  denominator  is  just  Equation  2)  evaluated  for  a  baseline 


29 


length  of  zero;  in  other  words,  the  antennas  are  adjacent.  This  shows 
that  the  correlation  coefficient  of  Brown  and  Twiss  is  just  the 
amplitude  of  the  fringe  visibility  function  derived  in  the  first 
section  of  this  chapter.  It  was  also  shown  at  that  time  that  the 
amplitude  of  the  fringes  generated  by  the  cross  correlation  function 
was  equal  to  the  fringe  amplitude  of  the  visibility  function. 

They  continue  in  their  paper  to  show  that  c  (£,t)  is  given  by 

F   2  (sj  +  F  .  2  (S  ) 
r  a   -n  -  cos   a    sm  v  a'   ,  .  ,    . 
c  l.*»cJ  ~  — ,  which  is  just  the  square  of  p  (£) . 

cos 

We  will  use  such  a  method  in  the  computer  to  determine  the  best  time 
shift  of  the  envelopes  of  bursts  received  at  two  different  stations  in 
order  to  determine  relative  arrival  times.  The  integration  period  used 
in  the  averages  of  the  functions  is  long  compared  with  any  time  variation 
in  the  arguments  of  the  voltages  sent  to  the  square-law  detectors. 
Figure  II. 3b  is  a  sample  record  of  data  taken  as  a  source  crosses  the 
antenna  pattern  of  the  individual  antennas. 

II-3.  The  Formation  of  an  Interferometer  Using  the 
Cross  Correlation  Function 

A  phase  interferometer  can  be  formed  using  data  that  is  tape 
recorded  at  the  antenna  locations.  The  receivers  must  heterodyne  the 
signal  from  the  observing  frequency,  in  this  case  18  MHz,  to  a  frequency 
compatible  with  the  tape  recorder.  The  data  must  be  digitized  and  is 
then  used  in  the  equation  for  the  cross  correlation  function 
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r  ,c  .  „,       E0  ct)  b  »(t  +  x)  i 

*   (.sx,C,tJ  =  _^__ u  .  The  form  actually  employed 

[E0Z  (t)]%  [E0Z  Ct  +  t)]* 


must,  of  course,  perform  the  integrations  discretely.  The  cross 

correlation  function   takes  the  form 

N 

£     E        (t)   E»       (t  +  t) 
C   Cs   ,C,t)   =     1=1     UJ °J .  The  variable  t,  is  a 


-N        N      -i% 

I      (E   2)   E   CE'.2) 
j=l   °3   j=l   °J  _ 

time  correction  factor  to  account  for  any  relative  timing  errors 
generated  by  such  things  as  transmission  time  for  timing  signals,  timing 
resolution,  etc.  The  method  involves  the  calculation  of  C  (s  ,£,t)  for 

A 

many  different  values  of  t,  and  the  selection  of  the  value  of  t  that 
maximizes  the  amplitude  of  C  (sx,5,t) .   Further  description  of  the 
procedure  follows  in  the  next  chapter. 

A  post- detection- correlation  interferometer  as  introduced  by  Brown 
and  Twiss  is  formed  by  first  square-law  detecting  the  heterodyned  data 
on  the  magnetic  tape.  After  smoothing,  the  envelopes  are  cross 
correlated.  A  value  for  the  best  relative  shift  of  the  envelopes  is 
determined  by  finding  that  value  of  t  that  maximizes  the  cross  correlation 
function    for  a  given  region  of  integration.  A  comparison  of  the  best 
shifts  found  for  the  two  methods  gives  a  means  to  detect  the  presence  of 
a  sweeping  beam  of  Jupiter  decametric  radiation. 

Occasionally  a  combination  of  circumstances  can  occur  in  which  a 
segment  of  data  that  has  been  interfered  is  so  short  in  time  that  the 
fringe  does  not  reach  maximum  value.  This  makes  the  determination  of  the 
fringe  amplitude  difficult  if  not  impossible.   If  the  local  oscillator 
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initial  phase  term  had  not  been  ignored  in  the  development  of  the  cross 
correlation   function,  C  (s  ,£,t),  a  function  of  the  form 

A 

C  (s  >5,n,t)  =  V  (s  ,£,n,t),  where  ri  is  the  initial  phase,  would  have  resulted. 

Equation  1)  in  section  II-l  would  have  the  form  V  (s. ,5,n,t)  =  V  (s  ) 

A  OX 

cos   [2tt  s      (<J>0  -  A(()0)   +  n  +  yt].     If  the  receiver  at  one  station  were 
arranged  such  that  two  equal  local  oscillator  frequencies  with  phases 
shifted  90°  were  injected  into  two  mixers,   the  signals  from  another 
station  cross  correlated  with  the  first  station  would  exhibit  fringes 
with  a  90°  phase  shift.     The  fringe  amplitude  V     (s  )   could  be  obtained 
by  taking  the  square  root  of  the  sum  of  the  squares  of  the  cross 
correlation     function       for  each  pair  of  signals.     The  derivation  is 
given  below: 

Vlst  pair  =  V0    ^   cos    V*  SX   C*0  "  AV   +  V  Tt] 


V2nd  pair  =  V0   <sx>   cos   •>  SX   C*0  "  AM  +  \  +  ^°  *  yt] 
V0    ^   Sin    [2TF  S^    C<,)n   "  A*n)    +  n*  +  Yt] 


[Vlst  pair  +  1„d  pai/  =  V0    <V    l™5*  {2«  SX   %  '  AV   +  "o  + 
yt}  +  sin2  {2tt  sx   (<j>0  -  A<j.Q)   +  n0  +  Yt}]3* 

therefore 

v0  (sj  =  [v2         .     CnJ  +  v2  Cnn  +  90°)]%  . 

0   X      1st  pair   0     2nd  pair   0 

For  a  strip  source  with  a  Gaussian  spatial  distribution  of 
radiators,  Equation  2) can  be  used  to  calculate  a  set  of  values  for  the 
amplitude  of  the  visibility  function,  V  (s  ) ,  if  the  source  is  symmetrical, 
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Since  it  has  been  shown  that  the  cross  correlation  function  and  the 
fringe  visibility  function  are  equal,  the  experimental  values  of  the 
amplitude  of  C(s  ,£,n,t)  should  lie  on  the  curve  representing  the  source. 
Figure  II. 4  is  a  graph  of  V  (sx)  plotted  versus  baseline  length  for 
sources  with  angular  widths  of  071,  0725,  and  075  of  arc.  Values  for 
the  amplitude  of  the  cross  correlation  function  determined  in  later 
chapters  will  be  compared  with  this  graph  to  determine  the  size  of  the 
source.  This  graph  is  adapted  from  Dulk  (1970). 
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CHAPTER  III 
THE  TAPE  RECORDING  INTERFEROMETER 

III-l.  The  Signal  Receiving  Systems 

In  1969  receiving  systems  as  diagrammed  in  Figure  III .1  were 
located  at  Western  Kentucky  University,  University  of  Florida  Radio 
Observatory,  Florida  Presbyterian  College,  and  the  Observatorio 
Radioastronomico  de  Maipu  in  Chile.  At  each  station  a  steerable, 
linearly  polarized  Yagi  antenna  provided  the  Jupiter  signal  to  a 
Collins  75S-1  communications  receiver.  The  receiver  used  crystal- 
controlled  local  oscillators  throughout  and  employed  a  mechanical 
filter  to  establish  the  overall  bandpass  of  the  receiver.  The 
mechanical  filter  had  a  455  kHz  center  frequency  and  a  2.1  kHz  band- 
width. A    50-ohm  resistor  serving  as  a  dummy  load  could  be  switched 
to  the  receiver  input  in  place  of  the  antenna.  A  WWV  time  standard 
receiver  was  operated  concurrently  with  the  system  to  provide  real 
time  synchronization  of  the  composite  time  channel.  A  data  channel 
multiplexer  received  the  audio  signals  from  the  Collins  and  WWV 
receivers  and  placed  them  on  one  channel  of  a  Magnecord  stereo  tape 
recorder.  The  multiplexer  placed  the  Jupiter  signal  on  the  tape  most 
of  the  time.  WW  was  keyed  in  for  two  seconds  every  minute  for  the 
double  tick  and  for  20  seconds  every  five  minutes  for  the  voice  announcement, 
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Figure  III.l  --  Diagram  of  1969  Receiving  System. 
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Before  the  voice  announcement  the  dummy  load  was  keyed  into  the  receiver 
input  in  place  of  the  antenna  for  a  few  seconds  at    5-minute  intervals 
to  serve  as  a  system  check.  A  crystal  oscillator  was  used  as  the  source 
of  a  stable  frequency  to  drive  a  divider  chain  to  provide  the  composite 
time  signal.  All  stations  had  the  one  pulse-per-second  and    60  pulses- 
per-second  outputs.  All  but  Western  Kentucky  University  used  the  960 
pulses-per-second  output.  A  600  pulses-per- second  output  was  used  by 
WKU.  The  composite  timing  signal  was  formed  by  pulse-height  coding 
the  tick.  The  one-second  pulse  was  tallest,  followed  by  the  l/60th- 
second  pulse  with  the  l/960th-second  (or  l/600th)  pulse  the  shortest. 
The  pulses  could  be  identified  on  a  chart  recording  or  oscilloscope  by 
this  means.  The  composite  time  signal  was  recorded  on  the  other 
channel  of  the  Magnecord  tape  recorder.  Real  time  synchronization  of 
the  composite  time  signal  was  accomplished  by  comparing  the  time  marks 
in  one  channel  with  the  onset  of  the  WWV  tone  or  double  ticks  that 
were  keyed  into  the  other  channel. 

The  processing  of  data  from  this  system  revealed  two  shortcomings. 
The  local  oscillators,  although  they  were  crystal  controlled,  were  not 
stable  enough.  One  could  not  predict  the  fringe  rate  generated  by  the 
local  oscillator  offset,  y,   since  this  offset  was,  in  itself,  not 
predictable.  The  other  shortcoming  occurred  when  two  stations  were 
interfered  at  a  time  when  their  fringe  was  passing  through  zero,  or  was 
a  value  other  than  a  maximum.   If  the  fringe  could  not  be  reconstructed 
for  a  long  enough  time  to  show  the  peak,  then  the  amplitude  of  the 
fringe  could  not  be  determined.  The  data  processed  for  this  disserta- 
tion from  1969  fortunately  had  a  duration  longer  than  one  fringe  cycle. 
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A  receiving  system  was  designed  for  the  1970  Jupiter  season  to 
overcome  these  shortcomings.  Figure  III. 2  is  a  block  diagram  of  the 
system  operated  at  the  University  of  Florida  Radio  Observatory.  A 
very  stable  crystal  in  a  double  oven  was  used  as  the  source  of  all 
frequencies  employed  in  the  system.  The  composite  time  signal  was 
produced  as  in  the  1969  system,  except  the  frequencies  were  changed  to 
1  pulse  per  second,  50  pulses  per  second  and  1000  pulses  per  second. 
The  time  signal  was  recorded  continuously  on  one  channel  of  a  four- 
channel  Viking  tape  recorder.  The  1.8  MHz  crystal  standard  also  drove 
a  frequency  synthesizer  that  provided  the  local  oscillator  injections 
to  the  receiver.  The  first  and  second  local  oscillator  frequencies 
were  produced  by  a  combination  of  multiplying  and  dividing.  The  second 
local  oscillator  frequency  was  switch  selectable  between  3.15  MHz  and 
4.05  MHz  to  provide  a  means  of  tuning  the  receiver  off  of  an  inter- 
fering station.  The  product  detector  local  oscillator  input  was 
produced  by  dividing  and  phase  shifting.  Two  product  detectors  were 
sent  the  signal  from  the  receiver  I.F.  stage.  Their  local  oscillator 
inputs  were  the  same  frequency,  450  kHz,  but  one  was  phase  shifted  90° 
with  respect  to  the  other.  The  interference  of  both  of  the  product 
detector  outputs  with  another  station  allows  the  reconstruction  of  the 
fringe  amplitude  without  having  either  combination  actually  reach  a 
fringe  peak.  The  receiver  was  designed  around  a  Collins  mechanical  filter 
of  the  type  used  in  the  earlier  system.  Two  different  bandwidths  were 
used,  2.1  kHz  and  6  kHz,  at  the  center  frequency  of  455  kHz.  The  R.F. 
preamp  and  mixers  were  made  using  dual-gate,  metal-oxide-silicon-field- 
effect  transistors  (MOSFET) .  They  provide  a  much  wider  dynamic  range 
and  superior  resistance  to  cross  modulation  due  to  adjacent  channel 
interference.  A  frequency  conversion  scheme  for  two  different  I.F. 
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Figure  III. 2  --  Diagram  of  1970  Receiving  System. 
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bandwidths  is  shown  in  Figure  III. 3.  A  detailed  block  diagram  of  the 
system  is  shown  in  Figure  III. 4.  The  receiver  located  at  the  other 
station  in  the  network  (Maipu,  Chile)  had  a  single  product  detector; 
otherwise  the  system  was  identical . 

The  real  time  synchronization  for  stations  in  the  United  States 
used  signals  from  WWV  and  LORAN-C.  The  data  multiplexer  keyed  WWV 
into  the  data  channel  every  minute  for  the  double  tick  and  every  five 
minutes  for  the  voice  announcement.  This  was  followed  by  LORAN-C  for 
two  seconds  to  provide  as  accurate  timing  as  possible.  The  VLF  trans- 
mission time  from  the  station  is  very  stable  and  predictable  for  ground 
wave.  Information  concerning  the  transmission  time  delay  for  WWV  time 
signals  came  from  Application  Note  52  by  the  Hewlett-Packard  Company 
(1965).   Information  pertaining  to  the  VLF  network  of  LORAN-C  came  from 
papers  by  L.  D.  Shapiro  (1968a  and  1968b).  The  inputs  to  the  four- 
channel  Viking  tape  recorder  consist  of:  channel  1  -  composite  time 
signal;  channel  2  -  continuous  WWV;  channel  3  -  multiplexed  data  and 
time  signals;  channel  4  -  data. 

Since  it  is  difficult  to  receive  any  of  the  LORAN-C  networks  in 
Chile,  they  were  provided  with  a  rubidium  frequency  standard  on  loan 
from  NRAO.  The  1  MHz  output  was  divided  down  to  1  pulse  per  second  and 
compared  with  NBS  time  at  the  NASA  satellite  tracking  station  near 
Santiago,  Chile.  The  atomic  standard  was  carried  while  still  running 
to  the  Jupiter  observatory  at  Maipu,  Chile,  to  calibrate  the  local 
standard  and  to  be  keyed  into  the  data  channel  every  minute  instead 
of  the  WWV  double  tick.  The  WWV  voice  announcement  was  keyed  in  every 
five  minutes  as  at  the  University  of  Florida  station. 
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The  antenna  system  at  Chile  consisted  of  a  polarimeter  that  was 
used  near  the  meridian  and  a  steerable  Yagi  for  sources  away  from  the 
meridian.  The  1970  data  processed  here  uses  the  linearly  polarized 
Yagi.  The  Florida  antenna  was  a  manually  steerable  polarimeter  or, 
on  the  night  when  the  data  was  recorded,  a  linearly  polarized  Yagi. 

The  network  has  since  been  enlarged  to  include  stations  at  Western 
Kentucky  University  (a  single  product  detector  model)  and  San  Antonio, 
Texas.  All  recording  in  the  network  in  1971  was  done  using  right  hand 
circular  polarization  except  at  WKU  (linear) . 

II 1-2.  The  Treatment  of  the  Data  on  the  Computer 

Interference  fringes  are  formed  for  various  pairs  of  stations  by 
calculating  the  cross  correlation  function   for  a  given  integration 
period  as  a  function  of  time.  In  order  to  use  the  computer  to  aid  in 
the  calculations,  it  was  necessary  to  present  the  data  to  it  in  a 
digital  form. 

The  data  tapes  were  digitized  using  an  analog-to-digital  converter 
as  a  data  source  for  an  IBM  1800  computer.  Two  forms  of  the  digitiza- 
tion program  were  used.  The  1969  data  was  digitized  by  passing  a  burst, 
that  had  been  located  earlier,  to  the  digitizer  using  the  original 
magnetic  tape.  The  operator,  listening  with  ear  phones  to  the  time 
channel,  counted  from  a  known  starting  time  (usually  the  WWV  minute  mark) 
to  the  second  tick  before  the  second  tick  that  preceded  the  part  to  be 
digitized.  He  pushed  the  START  button  on  the  computer,  which  searched 
for  the  next  second  mark.  Upon  finding  this  mark,  it  waited  for  a 
length  of  time  given  to  it  on  a  data  card.  After  the  wait,  the  computer 
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began  digitizing  until  it  filled  a  data  array  which  was  allotted  3520 
words.  The  time  channel  was  digitized  by  alternate  samples  with  the 
data  channel  and  stored  in  a  similar  array.  Both  arrays  were  transferred 
to  a  disk  and  plotted  on  an  incremental  plotter.  Figure  III. 5  shows 
a  composite  plot  of  the  digitized  data  and  time  channels  and  an 
expanded  plot  of  the  data  channel  for  Chile.  The  time  channel  is 
digitized  to  allow  the  exact  time  interval  of  digitization  to  be 
determined.  By  careful  comparison  with  the  chart  recordings  of  the 
raw  data  the  exact  starting  and  ending  times  of  the  digitization 
period  can  be  determined. 

The  data  from  1970  was  handled  a  bit  differently.  The  Computer 
Sciences  Division  of  the  University  of  Florida  Medical  Center,  where 
the  previous  work  had  been  done,  had  produced  a  general  purpose 
digitization  program  that  would  digitize  up  to  509,440  words  of  data 
at  speeds  in  excess  of  7  000  data  points  per  second  per  channel.  After 
the  data  was  discovered,  from  the  search  of  chart  records,  it  was 
slowed  down  by  a  factor  of  eight-to-one  and  presented  to  the  machine. 
Approximately  six  seconds  of  data  were  digitized,  with  points  being 
taken  alternately  from  the  data  channel  and  time  channel.  The  tape  had 
been  started  by  the  operator  and  digitization  begun  immediately.  A 
program  called  LOOK  was  provided  that  searched  the  data  for  the 
addresses  of  the  bursts  and  the  second  ticks.  The  interesting  bursts, 
which  were  previously  seen  on  chart  records,  were  then  located  and 
transferred  with  a  program  called  TRNFR  to  another  disk  for  storage 
and  plotting  on  the  incremental  plotter.  The  method  used  was  theoretically 
easier  than  that  used  the  previous  year,  but  tape  recorder  instabilities 
complicated  the  efforts. 
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Data  was  transferred  from  the  IBM  1800  on  cards.  The  remainder 
of  the  data  processing  was  done  on  an  IBM  360/65.  The  1970  data 
was  transferred  in  a  column  binary  format  in  order  to  require  the 
punching  of  fewer  cards  by  the  IBM  1800.  An  Assembler  program  was 
written  to  convert  the  column  binary  information  on  the  cards  to 
files  on  a  tape  that  could  be  read  by  a  FORTRAN  program  such  as 
CROCO,  which  is  described  later. 

In  order  for  interference  to  be  achieved  properly,  the  two  tapes 
from  a  pair  of  stations  must  present  data  to  the  multiplier  at  the 
same  rate  in  real  time.  Tape  stretch  and  other  inconsistencies  made 
some  correction  in  the  digitized  data  necessary.  A  program  called 
BADE  was  written  by  the  author  to  adjust  the  number  of  data  points 
in  a  table  to  represent  the  same  real  time  period  and  to  move  the 
zero  value  of  the  data  to  the  average  value  for  the  table.  The 
program  also  adjusted  each  data  table  to  begin  at  the  same  point  in 
real  time  by  correcting  for  time  signal  transmission  time  and  errors  in 
the  starting  time  of  digitization.   In  1969  a  data  table  of  2240  data 
points  was  prepared  by  BADE  such  that  the  first  point  in  each  table 
occurred  at  the  same  time.  This  left  an  error  due  to  the 
geometrical  path  length  differences  for  the  Jupiter  signal  and  some 
residual  timing  errors,  since  WWV  was  the  prime  source  of  real  time 
that  year. 

Two  programs  were  written  by  the  author  to  prepare  data  for  the 
intensity  interferometer  mode  and  the  envelope  correlation  experiment. 
The  program,  BASQ,  square- law  detected  the  data  in  each  of  the  data 
tables  prepared  by  BADE  and  integrated  it  over  20  data  points.  The  inte- 
gration period  was  sufficient  to  remove  the  audio  and  leave  the  slowly 
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time-varying  envelope.  A  post-detection-correlation  interferometer 
of  the  Brown  and  Twiss  type  was  formed  with  this  data.  The  other 
program,  BABS,  took  the  absolute  value  of  the  BADE  data,  integrated 
over  20  data  points,  and  adjusted  the  baseline  of  the  data  to  the 
average.  It  played  the  part  of  a  full-wave,  linear  detector  with 
integrator,  followed  by  a  coupling  capacitor.  BABS  provided  the 
data  to  examine  the  relative  arrival  times  of  the  fine  structure 
in  an  L-Burst. 

The  program,  called  CROCO,  was  the  main  analysis  tool  used  on 
the  data.  Its  purpose  was  to  calculate  the  cross  correlation  function, 
C(sx,5,n,t) .  It  was  written  by  the  author  to  be  very  flexible  in  its 
choice  of  integration  times,  number  of  values  of  C  calculated  for  a 
data  table,  and  where  (in  a  data  table)  it  started  and  ended  calcula- 
tions of  C.  CROCO  performs  the  calculation  of  the  cross  correlation 
function, 
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in  two  parts.  The  numerator  is  formed  by  summing  the  product  E 

E0  -+ITAU  over  N  vaIues  of  j ,  where  N  is  the  selected  integration  time 

and  ITAU  is  the  initial  shift  of  the  two  data  tables  with  respect  to 

each  other.  The  denominator  is  formed  by  computing  the  sum  of  the 

square  of  E   .  and  the  sum  of  the  square  of  E'       .  The  two  sums 
0  3  0  j+ITAU 

are  multiplied  together  and  the  square  root  is  taken.  The  cross 
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correlation  function   is  formed  and  given  the  name  PXC(J),  an 
element  of  an  array.  The  denominator  is  saved  and  stored  as  PNM. 
The  program  steps  LDELT  data  points  further  into  the  data  tables  and 
calculates  another  value  for  C  and  the  denominator.  After  calculating 
three  values  of  C,  a  line  of  print  is  generated.  The  line  of  print 
consists  of  the  denominator  (PNM),  three  values  of  C  (XC0R1,  XC0R2, 
XC0R3)  and  the  line  number.  In  addition,  the  program  assigns  an 
asterisk  to  one  of  the  100  remaining  spaces  in  the  line  at  a  position 
that  is  proportional  to  the  value  of  the  first  C  calculated.  Since  C 
ranges  from  -1  to  +1,  each  space  represents  a  change  in  C  of  .02. 
The  process  is  repeated  until  the  index  for  the  selection  of  the  data 
points  to  be  multiplied  reaches  NDPS,  which  is  the  number  of  data 
points  in  each  table.  The  program  increments  the  relative  shift 
variable  K  by  KDELT  and  repeats  the  calculation.  When  K  has  the 
value  JTAU,  the  program  will  repeat  the  cycle  once  more  and  stop. 

Some  essential  quantities  which  must  be  provided  by  the  user  are 
the  following: 

1.  IDATE(l),  IDATE(2),  AND  IDATE(3)  are  the  numbers  of  the  month, 
day  and  year,  respectively,  e.g.,  1,  2,  69,  on  which  the  burst  occurred, 

2.  IFIL1  and  IFIL2  are  used  to  denote  the  portion  of  an  input 
array  in  which  data  from  each  station  has  been  placed.  The  initial 
array  usually  must  be  changed  for  different  batches  of  data. 

3.  JHOUR  and  XSEC  are  the  hour  and  second  at  which  the  burst  of 
interest  occurred. 

4.  ITAU  and  JTAU  are  the  beginning  and  ending  relative  shifts. 
They  are  calculated  by  the  formulas  ITAU  =  (initial  time  shift)  x 
(digitization  rate)  and  JTAU  =  (final  time  shift)  x  (digitization  rate) . 
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If  the  two  shifts  are  on  the  same  side  of  zero-relative-shift,  ITAU 
is  given  a  negative  sign. 

5.  N  is  the  integration  time  in  data  points.  N  =  (integration 
time)  x  (digitization  rate) . 

6.  KDELT  is  the  increment  of  shifting  in  data  points. 

7.  LDELT  is  the  increment  of  stepping  through  the  data  tables. 

8.  DPS  is  the  digitization  rate  in  data  points  per  second. 

9.  NDPS  is  the  number  of  data  points  in  each  data  table.  All 
tables  are  assumed  to  contain  the  same  number  of  data  points. 

The  program  is  run  on  data  from  a  pair  of  stations  for  a  given 
range  of  relative  shifts.  The  best  shift  is  selected  by  noting  the 
shift  that  makes  the  fringes  have  the  greatest  amplitude  and  have  the 
smoothest  curve. 

A  modification  of  CROCO  was  made  that  incorporated  statements 
to  control  an  incremental  plotter.  CROCOCAL  produced  the  graphs  of  the 
fringes  shown  in  a  number  of  the  figures.  It  was  hoped  that  it  would 
be  inexpensive  enough  to  use  as  an  analytic  tool  in  place  of  CROCO, 
since  the  smoothness  of  the  curves  is  much  easier  to  see.  Unfortunately, 
it  turned  out  to  be  anything  but  inexpensive.  A  listing  of  a  combina- 
tion version  of  BADE  and  BASQ  is  included  in  the  Appendix, 
followed  by  a  listing  of  CROCO. 

A  refinement  in  the  method  of  setting  the  digitization  rate  was 
added  during  the  processing  of  the  1970  data.  An  integrated  circuit 
phase-lock  detector  made  by  Signetics  (NE  565)  was  used  to  lock  an 
oscillator  to  the  15  KHz  tone  that  had  been  placed  on  the  time  channel. 
Since  the  tone  had  been  derived  from  the  crystal-frequency  standard, 
it  remained  in  phase  with  the  timing  marks  on  this  channel.  The  analog- 
to-digital  converter  in  the  IBM  1800  Data  Acquisition  System  was  driven 
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from  the  synchronized  oscillator.  The  synchronization  of  the  converter 
eliminates  any  change  in  the  effective  digitization  rate  due  to 
changes  in  tape  speed  (caused  either  by  differences  in  tape  recorder 
speeds  or  wow  and  flutter  in  a  given  tape  recorder) .  The  time  channel 
of  the  digitized  data  was  printed  by  the  IBM  1800  in  such  a  form  as 
to  make  detection  of  any  loss  of  coherence  between  the  timing  channel 
and  the  phase-locked  oscillator  obvious.  For  further  checking  the 
data  channel  was  plotted  on  the  incremental  plotter.  The  data  was 
punched  on  cards  in  a  column  binary  format  for  transportation  to  the 
IBM  360  for  use  in  CROCO. 


CHAPTER  IV 
AN  ANALYSIS  OF  THE  DATA  RECEIVED  IN  1969 

IV- 1.  Spatial  Orientation  of  the  Receiving  System 

The  Jovian  decametric  burst  that  is  analyzed  in  this  chapter  was 
received  on  the  morning  of  2  January  1969  at  08h22m14s  UT  by  the  four 
stations  that  were  members  of  the  interferometer  network  that  year. 
The  stations  were  located  at  Bowling  Green,  Kentucky  (lat.  36°57'  N, 
long.  86°25*  W) ,  Old  Town,  Florida  (lat.  29°31'50"  N,  long,  82°1' 
55"  W),  St.  Petersburg,  Florida  (lat.  27°45l45"  N,  long.  82°38»  W) , 
and  Maipu,  Chile  (lat.  33°23'50"  S,  long.  60°42'12"  W) .  The  primary 
timing  source  was  the  National  Bureau  of  Standards  radio  station  WWV, 
located  in  Boulder,  Colorado  (lat.  40°40'48"  N,  long.  105°2»25"  W) , 
as  shown  in  Figure  IV. 1.  The  local  one-second  ticks  were  delayed  by 
approximately  one  second  from  the  one- second  tick  transmitted  by  WW 
due  to  the  transmission  time  delay  and  the  manual  method  used  for 
aligning  the  local  tick  with  that  from  WW .  The  total  delays  for  the 
stations  were  1.0016  seconds  for  Bowling  Green,  0.4202  seconds  for 
Old  Town,  0.8917  seconds  for  St.  Petersburg,  and  0.2555  seconds  for 
Maipu.  The  actual  transmission  delay  from  WWV  ranged  from  6.1  milli- 
seconds for  Bowling  Green  to  about  33  milliseconds  for  Maipu.  The  latter 
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Figure  IV. 1  --  Western  Hemisphere  Map  Showing  the  Locations  of  the 
Receiving  Stations  and  WWV.  MAIPU  Is  near  Santiago, 
Chile. 
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time  is  uncertain  by  at  least  2  milliseconds  since  it  is  not  known,  a 
£riori,  how  many  hops  were  required  for  transmitting  the  timing  signal 
that  was  received  at  Maipu. 

Figure  IV. 2a  shows  the  arrangement  of  the  earth,  Jupiter,  and 
the  sun  on  2  January  1969  as  seen  from  the  North  Pole  of  the  ecliptic. 
The  right  ascension  of  Jupiter  was  12n22m19?2  and  its  declination  was 
-0°59'28".  The  right  ascension  of  the  sun  on  that  date  was  approximately 
18  41  or  280°  east  of  the  Vernal  Equinox.  The  North  Pole  of  the  earth 
is  pointing  to  the  right  as  the  observer  looks  along  the  line-of -sight 
from  the  earth  to  Jupiter,  since  this  is  the  northern  hemisphere  winter. 

The  receiving  stations  will  be  referred  to  by  the  initials  of  the 
university  radio  observatory  that  operated  them.  They  were  Western 
Kentucky  University  (WKU)  at  Bowling  Green,  University  of  Florida  Radio 
Observatory  (UFRO)  at  Old  Town,  Florida  Presbyterian  College  (FPC)  at 
St.  Petersburg,  and  Observatorio  Radioastronomico  de  Maipu  of  the 
Universidad  de  Chile  (MAIPU)  at  Maipu,  Chile.  Table  IV. 1  gives  the 
baseline  lengths  of  pairs  of  these  stations  for  the  baselines  projected 


Station  Pair 
UFRO  -  FPC 
WKU  -  UFRO 
WKU  -  FPC 
FPC  -  MAIPU 
UFRO  -  MAIPU 
WKU  -  MAIPU 


Table  IV. 1 

DEW   (km) 

DNS    (km) 

D    (km) 

83.1 

166 

185.8 

44.7 

697 

698 

89.6 

896 

905 

2000 

6520 

6820 

1913 

6690 

6980 

1875 

7380 

7620 

Vernal  Equinox 
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Figure  IV. 2  --  (a)  View  of  Earth,  Jupiter  and  the  Sun  from  Above  the 

North  Pole  on  2  January  1969,  0>)  D   and  D^  Projected 

NS      BW 
on  the  Plane  Perpendicular  to  the  Line-of-Sight  to 

Jupiter. 
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on  a  plane  perpendicular  to  the  line-of-sight  to  Jupiter  from  the  earth. 
The  variable  DEW  is  the  length  of  the  baseline  projected  on  a  line  that 
is  parallel  to  the  equatorial  plane  and  is  perpendicular  to  the  line- 
of-sight.  The  variable  DNS  is  the  length  of  the  baseline  projected  on 
a  line  parallel  to  the  north-south  axis  of  the  earth  and  perpendicular 
to  the  line-of-sight  to  Jupiter.  The  variable  D  is  the  total  baseline 
projected  on  a  plane  perpendicular  to  the  line-of-sight  to  Jupiter. 
The  effective  baseline  lengths,  measured  in  wavelengths  of  the  18  MHz 
received  signal  of  the  various  interferometers  formed  by  all  pairs  of 
stations,  range  from  11,150  to  456,000. 

For  the  experiment  involving  the  effect  of  the  solar  wind  on  the 
envelope  of  the  L-Burst  (as  proposed  by  Douglas  and  Smith) ,  it  is 
necessary  to  know  the  effective  length  of  the  baseline  as  projected  on 
the  plane  of  the  ecliptic.   It  is  the  interaction  of  the  solar  wind 
with  the  radiation  travelling  from  Jupiter  along  the  line-of-sight, 
which  lies  in  the  ecliptic  plane,  that  is  proposed  to  be  the  cause  of 
the  shape  of  the  L-Burst  envelope.  Figure  IV. 2b  shows  the  appearance 
of  the  WKU  -  MAIPU  baseline  components  D   and  D   as  seen  by  an 
observer  at  the  center  of  the  earth  looking  toward  Jupiter.  The  direction 
sense  for  east  and  west  is  that  used  by  an  observer  at  the  location  of 
the  stations,  i.e.,  on  the  side  of  the  earth  toward  Jupiter  in  Figure 
IV. 2a.  The  angle,  a,  is  the  projection  of  the  polar  angle  of  the  earth 
(23%°) -on  the  plane  perpendicular  to  the  line-of-sight  to  Jupiter.  Its 

value  is  given  by  tan  a  =  tan  23V  cos  (RA  -180°),  where  RA  is  the 

J  J 

right  ascension  of  Jupiter  in  degrees.  On  2  January  1969  the  right 
ascension  of  Jupiter  was  185?5.  The  value  of  a  was  then  23°4.  The 
east-west  baseline,  BEW,  as  projected  on  the  plane  of  the  ecliptic, 
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is  given  by  B   =  DNg  sin  a  +  D   cos  a  for  the  orientation  shown 
in  Figure  IV. 2b. 

IV-2.  The  Nature  of  the  Received  Signal 

According  to  a  classification  of  Jovian  noise  storm  activity- 
given  by  Douglas  and  Smith  (1967)  and  attributed  to  Gallet,  a  burst 
with  a  duration  of  from  0.1  second  to  5  seconds  is  an  L-Burst. 
Inspection  of  Figures  IV. 3  and  IV. 11  will  reveal  that  the  burst  under 
study  here  is  about  0.130  seconds  long;  hence,  by  strict  application 
of  the  classification  scheme,  it  is  technically  an  L-Burst.  Recall 
that  the  largest  time  marks  that  are  visible  in  the  figures  occur  at 
l/60th  second  intervals.  This  burst  will  henceforth  be  referred  to 
as  an  L-Burst,  even  though  an  added  requirement  for  a  burst  to  be  so 
classified  will  be  shown  as  a  result  of  the  analysis  in  this  chapter. 

As  can  be  seen  in  Figure  IV. 3,  the  signals,  as  received  at  the 
several  stations,  were  not  all  of  the  same  strength.  An  important 
requirement  for  the  use  of  the  signals  in  the  interferometer  is  their 
signal-to-noise  ratio.  This  number  is  gotten  by  comparing  the  voltage 
ratio  of  the  peaks  of  the  burst  with  the  peaks  of  the  noise  outside  of 
the  burst.  Expressed  in  decibels,  the  (S+N)/N  ratios  are  the  following: 
WKU  -  14  db,  FPC  -  9.5  db,  MAIPU  -  9  db,  and  UFRO  -  6  db.  The  signal- 
to-noise  ratio  for  UFRO  is  too  low  to  make  the  fringes  detectable. 
Another  factor  in  the  form  of  the  relative  frequency  effect  of  the 
local  oscillator  at  UFRO  made  the  fringes  almost  indistinguishable 
from  other  noise  effects.  This  required  the  exclusion  of  the  UFRO  data 
from  the  calculations  of  the  cross  correlation  function. 
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Figure  IV. 3  --  Actual  Chart  Records  of  the  Heterodyned  Data  from  WKU  (a), 
FPC  (b),  UFRO  (c)  and  MAIPU  (d)  for  the  L-Burst  of 
2  January  1969.  Alternate  Lines  Are  the  Timing  Channels 
for  the  Corresponding  Stations. 
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IV-3.     Measurement  of  the  Size  of  the  Source  of 
Jovian  Decametric  Radiation 

It  was  shown  in  Chapter  II  that  the  brightness  distribution  of  a 
source  could  be  calculated  from  the  fringe  visibility  function  V  fs  1 

by  the  following  equation:     B(<fr  J   =  S     /     V.  (s  Je~l2lTSA(VAVds   . 

In  that  chapter  it  was  also  shown  that  the  visibility  function  was 
equal  to  the  amplitude  of  the  cross  correlation  function,  C(s    ,5,t). 

A 

Therefore,  in  order  to  find  the  angular  width  (i.e.,  twice  the  value 
of  <j>0  for  which  B C*QD  is  just  zero)  of  the  source,  one  must  calculate 
the  cross  correlation  function  for  the  best  relative  shift  of  the  data 
from  all  pairs  of  stations.  The  best  relative  shift  is  determined  by 
the  relative  shift,  t,  in  Equation  4)  in  Chapter  III  that  maximizes  the 
amplitude  of  the  curve  of  the  cross  correlation  function  plotted 
against  time.  The  reason  t  was  introduced  was  to  make  allowance  for 
the  relative  inaccuracy  of  the  knowledge  of  real  time  at  each  of  the 
stations  (due  to  real  time  being  received  from  WWV)  and  for  the 
difference  in  path  length  that  the  signal  from  Jupiter  had  to  travel 
to  the  various  stations. 

The  program  CROCO  was  used  on  the  heterodyned  data,  which  had 
previously  been  adjusted  in  time  and  length  by  the  program  BADE,  to 
calculate  and  make  a  plot  of  the  cross  correlation  function  as  a  function 
of  time.  The  plot  was  in  the  form  of  a  series  of  asterisks  on  the  computer 
printout,  as  is  described  in  Chapter  III.  After  the  best  shift  was 
discovered,  the  program  CROCOCAL  plotted  the  curves  that  are  shown  in 
Figures  IV. 4,  IV. 5,  and  IV. 6.  The  individual  figures  are  composites 
of  three  pages  of  plotting  of  the  output  of  CROCOCAL  done  on  a 
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Calcomp  plotter.  The  first  two  lines  are  a  plot  of  the  data  as 
received  by  the  stations  (the  heterodyned  signal) ,  hence  can  be 
compared  directly  with  Figure  IV. 3.  The  third  line  in  each  figure  is 
the  cross  correlation  function  plotted  against  time,  and  the  fourth 
line  is  the  denominator  of  C(s  ,5,t)  which  serves  as  a  normalizing 

A 

factor.  The  integration  time,  N,  is  2.1  milliseconds.  The  peaks  in 
the  normalization  curve  indicate  where  especially  strong  regions  of 
the  burst  are  located.  The  time  axis  is  calibrated  in  data  points, 
such  that  one  data  point  is  0.07  milliseconds.  Time  is  counted  from 
the  beginning  of  the  burst,  or  08h22m14§3  UT. 

Figure  IV. 4  shows  the  data  from  WKU  compared  with  that  from 
MAIPU.  Figure  IV. 5  compares  data  from  WKU  and  FPC,  while  Figure  IV. 6 
compares  data  from  FPC  and  MAIPU.  Table  IV. 2  gives  the  amplitude  of 
the  cross  correlation  function  near  the  two  selected  regions  of 
high  signal  level  and  shows  the  fringe  rate  in  this  vicinity.  As  was 
shown  in  Chapter  II,  the  fringe  rate  is  due  to  the  relative  motion 
of  the  source  with  respect  to  the  far  field  pattern  of  the  inter- 
ferometer, together  with  the  relative  local  oscillator  offset  frequency. 
The  fringe  rate  due  to  the  motion  of  the  source  through  the  inter- 
ferometer pattern  for  the  longest  east-west  baseline  is  on  the  order  of  0.1 
fringe/second.  Most  of  the  fringe  rate  in  Table  IV. 2  can  be  attributed 

Table  IV. 2 

Stations     |c|  near  P.P.  880  |c|  near  P.P.  1240      F.R. 

WKU  -  FPC       0.90  ±  14%  0.92  ±  14%          107  Hz 

FPC  -  MAIPU      0.80  ±20%  0.87  ±  20%          216  Hz 

WKU  -  MAIPU     0.95+16%  0.80  ±  16%          103  Hz 
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to  the  local  oscillator  offsets.  It  is  encouraging  to  note  that  the 
sum  of  two  of  the  fringe  rates  (the  first  and  the  third)  is  very 
nearly  equal  to  the  fringe  rate  of  the  other  pair  of  stations  (the 
second)  as  is  predicted  by  the  theory  of  operation  of  the  super- 
heterodyne receivers  used  in  this  experiment. 

It  is  clear  from  Table  IV. 2  that  the  amplitude  of  the  cross 
correlation  function   does  not  decrease  with  baseline  length.  The 
variation  in  |c|  is  closely  related  to  the  signal-to-noise  ratio  of 
the  signals  used  to  calculate  it.  This  effect  is  also  noticed  in  the 
relative  "fuzziness"  of  the  curve  for  C(sx,£,t)  in  regions  where  the 
signal  is  not  strong,  as  is  reflected  by  the  normalization  curve. 
The  immediate  conclusion  from  this  data  is  that  the  longest  inter- 
ferometer baseline  did  not  resolve  the  source  of  this  burst.  There  is 
another  interpretation  of  this  result,  however.   If  the  source  is 
coherent,  the  assumptions  made  concerning  the  random  nature  of  the 
signal  as  a  function  of  both  time  and  position  on  the  source  are 
invalid.  The  result  is  that  the  fringe  visibility  function,  VQ 
cannot  be  used  to  reconstitute  the  brightness  distribution.  This 
consequence  will  be  considered  in  the  concluding  portion  of  the 
dissertation.  At  this  point,  it  can  be  concluded  that  the  angular 
width  of  the  source  is  less  than  the  width  of  one  beam  lobe,  if  the 
source  is  incoherent.  This  angle,  for  the  longest  baseline,  is  given 
by  6  =  -  radians.  For  the  WKU  -  MAIPU  baseline,  D  is  7620  km, 
making  6  have  a  value  on  the  order  of  h   second  of  arc.  The  corresponding 
linear  distance  on  Jupiter  is  1880  km. 

A  smaller  upper  limit  on  the  size  of  the  source  can  be  found  by 
using  the  family  of  curves  for  the  amplitude  of  the  fringe  visibility 
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function,  VQ(sx),  that  ;is  given  in  Figure  II. 4.  Using  the  inter- 
pretation that  low  values  of  |c|  are  due  to  uncorrelated  noise 
contamination  when  the  signal  strength  is  weaker,  it  is  seen  that  the 
data  points  form  a  curve  that  lies  above  the  curve  for  a  071  of  arc 
Gaussian  strip  source.  This  implies  that,  if  the  source  is  incoherent, 
its  size  is  less  than  0V1  of  arc  or  a  linear  distance  of  490  km 
on  Jupiter. 

An  experiment  was  conducted  using  a  one-bit  correlation  scheme  to 
compare  the  method  with  that  already  given.  Weinreb  (1963)  and  Cooper 
(1970)  have  demonstrated  the  utility  of  various  one-  and  two-bit 
correlation  techniques  in  the  field  of  radio  astronomy.  The 
sensitivity  of  the  one-bit  technique  is  64%  that  of  the  continuous 
multiplying  correlator  used  above.  The  various  two-bit  correlation 
methods  have  a  sensitivity  on  the  order  of  85%  compared  with  the 
continuous  multiplying  correlator.  It  was  felt  that  by  converting  the 
data  to  the  form  in  which  plus  one  represents  excursions  of  the  data 
voltage  above  the  zero  axis  and  minus  one  those  below  the  axis,  any 
contribution  of  the  envelope  to  the  cross  correlation  function 
would  be  eliminated.   It  should  be  recalled  that  there  are  many 
physical  modifications  that  can  occur  to  the  envelope  of  the  Jovian 
burst,  e.g.,  solar  wind  interaction,  interaction  of  a  Faraday  rotated 
plane-polarized  wave  and  an  antenna  polarized  in  another  sense,  etc. 
Any  frequency  effect,  having  been  introduced  at  the  source,  is  useful 
for  the  purpose  of  giving  information  about  the  source  itself.  Figure 
IV. 7  is  a  composite  of  the  fringes  generated  by  interfering  signals 
that  have  been  treated  in  the  one-bit  method  from  the  stations  WKU 
and  MAIPU.  The  one-bit  representations  of  the  signals  from  WKU  and 
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MAIPU  are  shown  on  lines  one  and  two,  respectively.  The  third  line  is 
the  cross  correlation  function  plotted  against  time,  with  N  =  2.1 
milliseconds,  and  the  fourth  line  is  the  normalization.  The  normaliza- 
tion has  a  constant  value  since  the  amplitude  information  in  the 
original  signals  has  been  removed.  This  is  one  of  the  attractive 
features  of  the  method  from  an  instrumentation  point  of  view.  The 
hardware  for  the  correlator  is  relatively  easy  to  implement  with 
digital  logic  integrated  circuits.  The  time  delay,  t,  can  be  handled 
by  shift  registers,  for  example.  The  cross  correlation  function 
plotted  in  Figure  IV. 7  shows  the  fringes  that  were  seen  in 
Figure  IV. 4,  which  is  a  plot  of  the  same  data  without  the  one- 
bit  representation.  Some  regions,  corresponding  to  highly  intense 
radiation  in  Figure  IV. 4,  show  a  fringe  amplitude  on  the  order  of 
0.9.   It  should  be  noticed  in  Figure  IV. 7  that  the  fringe  amplitude 
decreases  away  from  the  center  of  the  burst.  This  artifact  is  caused 
by  the  time  scaling  having  some  residual  error  which  causes  the  data 
at  the  beginning  and  end  of  the  burst  not  to  have  exactly  the  same 
best  shift.  This  is  removed  in  Figure  IV. 4  by  plotting  the  first  and 
last  thirds  of  the  data  at  the  proper  best  shift  for  that  section  of 
the  data.  The  one-bit  correlator  is  seen  to  depend  very  sensitively 
on  relative  shift  of  the  data.  This  feature  will  be  used  in  the 
section  of  Chapter  IV  dealing  with  the  beam-sweep  experiment.  The 
principal  shortcoming  of  the  one-bit  correlator  is  that  it  must  integrate 
over  a  longer  period  of  time  to  overcome  the  loss  of  information  due 
to  the  one-bit  representation.   In  normal  radio  astronomical  appli- 
cations there  is  no  difficulty  encountered  with  this  shortcoming  since 
the  sources  that  are  studied  do  not  change  in  intensity  over  times 
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that  are  long  with  respect  to  a  typical  integration  period.  The 
Jovian  burst  envelope,  however,  changes  with  time  and  its  shape  is 
not  preserved  from  one  burst  to  another.  The  S-burst,  especially, 
requires  the  use  of  the  continuous  multiplying  correlator. 

In  Chapter  II,  it  was  shown  that  the  Brown  and  Twiss  post-detector- 
correlation  interferometer  could  be  formed  using  the  same  form  of  the 
cross  correlation  function  as  was  used  for  the  phase  interferometer. 
The  data  from  the  four  stations  was  full-wave  detected  in  the  computer 
by  the  program  BABS.  CROCO  was  run  using  this  data  in  the  same 
manner  as  before.  The  best  shift  for  each  pair  of  stations  was 
found,  and  the  cross  correlation  function  was  plotted  by  CROCOCAL. 
Figures  IV. 8,  IV. 9,  and  IV. 10  show  the  results  of  these  calculations 
for  the  three  pairs  of  stations  used  in  the  phase  interferometer. 
In  each  figure,  lines  one  and  two  are  the  detected  data  from  the 
stations.  The  full-wave  detected  data  was  filtered  for  an  integra- 
tion time  of  1.4  milliseconds.  The  integration  time,  N,  in  the 
cross  correlation  function  calculation  was  70  milliseconds,  which  is 
approximately  the  length  of  the  detail  structure  in  the  data.  The 
third  line  is  the  cross  correlation  function,  which  no  longer 
exhibits  fringes.  The  fourth  line  is  the  normalization.  Table  IV. 3 
summarizes  some  values  of  the  cross  correlation  function  near 


Stations 
WKU  -  FPC 
FPC  -  MAIPU 
WKU  -  MAIPU 


Table  IV. 

3 

C2  near  D.P. 

880 

C2 

near  D.P. 

1240 

0.84 

0.97 

0.64 

0.95 

0.90 

0.97 
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regions  of  especially  strong  signal  (the  same  regions  as  were  shown 
in  Table  IV. 2).  The  symbol  C2  means  that  this  cross  correlation 
function  is   to  be  interpreted  as  the  square  of  the  C  found  using 
the  phase  interferometer.  It  is  seen  that  the  values  obtained  by  the 
two  methods  are  essentially  the  same.  Since  the  cross  correlation 
function,  C  ,   also  does  not  decrease  with  increasing  baseline  length 
it  can  be  concluded  that,  if  the  source  is  incoherent,  it  is  not 
resolved.  The  beam  lobe  angle  is  the  same  as  for  the  phase  inter- 
ferometer. If  the  source  of  this  burst  was  incoherent,  it  must  have 
been  less  than  071  of  arc  in  angular  width. 

IV-4.   Burst  Arrival  Time  Experiment 

Figure  IV. 11  is  a  composite  of  three  chart  records  depicting 
the  envelope  of  the  burst  as  it  was  received  at  WKU,  FPC  and  MAIPU. 
The  envelopes  have  been  aligned  in  real  time,  i.e.,  local  tick 
misadjustment  and  IVWV  time  delay  have  been  removed.   It  is  clear  that 
this  burst  arrived  at  all  of  the  stations  at  essentially  the  same  time. 

A  much  more  accurate  technique  can  be  used  than  that  given  above. 
The  normal  method  to  measure  the  burst  arrival  time  consists  of 
providing  stable  local  time  marks  on  the  data  tape  that  have  been 
accurately  calibrated  in  real  time.  The  calibration  procedure  is  very 
difficult  and  has  only  recently  become  practicable  for  widely 
separated  stations.  The  method  usually  involves  flying  an  atomic 
oscillator  by  plane  to  the  various  stations  or  comparison  with  a 
continuously  calibrated  time  standard  in  an  earth  satellite.   In  the 
novel  technique  for  measuring  burst  arrival  time  given  here,  Jupiter's 
signals  themselves  are  used  to  calibrate  the  local  clocks. 
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Figure  IV. 11 


•-  Actual  Chart  Records  of  the  Detected  Data  from  WKU   (a) , 
FPC   (b)     and  MAIPU  (c)   for  the  L-Burst  of  2  January  1969, 
Records  Are  Laterally  Adjusted  Such  as  to  Be  Aligned 
in  Real  Time. 
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The  property  of  the  Jovian  signals  that  is  used  can  be  demon- 
strated by  considering  a  model  consisting  of  a  frequency  modulated 
transmitter  driving  a  directive  antenna  that  is  located  in  the 
vicinity  of  Jupiter  (it  is  only  necessary  that  the  station  be  a  long 
distance  away  and  in  the  line-of -sight) .  If  two  receiving  stations 
are  located  such  that  the  beam  of  the  antenna  can  sweep  across  them, 
the  output  of  the  receivers  will  be  as  shown  in  Figure  IV. 12.   It  is 
assumed  that  the  station  that  is  not  in  the  center  of  the  beam  can 
still  receive  some  signal.  The  receivers  are  superheterodynes  and 
have  product  detectors,  as  do  those  in  our  interferometer.  Figure 
IV. 12  has  been  drawn  in  such  a  way  that  the  phase  changes  in  the 
frequency  modulation  are  aligned  in  each  of  the  signals,  e.g.,  the 
feature  at  time  "7".  The  frequency  modulation  features  are  produced 
at  the  transmitter  and,  except  for  random  Doppler  shifting,  cannot  be 
altered  during  transmission  after  the  signal  leaves  the  transmitter. 
The  envelopes,  on  the  other  hand,  are  determined  by  the  location  of 
the  receiver  with  respect  to  the  antenna  beam  of  the  transmitter.   If 
the  beam  sweeps  from  one  station  to  the  other,  there  will  be  a 
noticeable  difference  in  the  time  relationship  of  the  frequency 
modulated  "carrier"  and  the  envelope  when  the  signals  received  at  each 
station  are  compared.  The  frequency  modulation  feature  at  time  "7" 
has  a  high  amplitude  at  station  two,  while  it  is  weaker  at  station 
one.  If  the  signal  underneath  the  envelope  from  one  station  is  cross 
correlated  with  that  from  the  other,  the  relative  shift  for  the  greatest 
amplitude  of  the  cross  correlation  function    will  be  zero  for  this 
example.  If  the  envelopes  are  cross  correlated,  it  is  clear  from 
Figure  IV. 12  that  the  envelope  of  station  two  would  have  to  be  shifted 


73 


• 


• 


/i 


«=1 

N. 


N 


\ 


M 


--LJI 


r\ 


L> 


f 


U 


I/' 

STATION    2 


I J I I L 


j- 


J L 


J L 


4  6 

TIME 


8 


10 


Figure  IV. 12  -• 


Plot  of  the  Amplitude  of  a  Constant  Intensity  Signal 
as  Received  at  Two  Stations  as  the  Beam  of  Radiation 
Sweeps  from  Station  1  to  Station  2.  Frequency 
Modulation  Features  Are  Visible  under  the  Envelopes. 
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forward  by  approximately  3  time  units  to  obtain  the  maximum  cross 
correlation  function. 

It  is  proposed  by  a  number  of  investigators,  including  Goldreich 
and  Lynden-Bell  (1969),  that  the  decametric  radiation  from  Jupiter 
is  highly  beamed.  A  burst  envelope  received  on  earth,  then,  may  be 
shaped  by  the  way  that  the  beam  sweeps  across  the  receiving  station, 
if  it  sweeps  rapidly  enough.  The  frequency  modulation  of  the  Jovian 
signal  is  due  to  the  source  being  a  producer  of  random  noise  of 
moderate  bandwidth.  It  is  the  interaction  of  this  noise  spectrum 
with  the  I.F.  filter  in  the  receiver  that  produces  the  wave  forms 
shown  in  Figure  IV. 3.  The  filters  in  all  of  the  receivers  in  the 
interferometer  are  identical,  hence  their  effect  on  the  same  noise 
spectrum  is  the  same.  There  is  only  one  best  shift  for  the  heterodyned 
data  due  to  the  random  nature  of  the  noise  spectrum  produced  by  the 
Jovian  decametric  source.  The  local  clocks  can  be  calibrated  with 
each  other  by  using  the  correction  produced  by  finding  the  best  shift 
for  each  pair  of  stations.  The  envelopes  of  the  signals  can  be 
cross  correlated  to  find  the  best  shift  for  them.  The  best  shifts  for 
the  heterodyned  data  and  the  envelopes  from  a  given  pair  of  stations 
can  be  compared  to  determine  the  rate  and  direction  of  motion  of  a 
sweeping  beam. 

Figure  IV. 13  is  a  plot  of  a  portion  of  the  burst  as  received  at 
WKU  and  MAIPU.  The  signal  is  especially  strong  in  this  part  of  the 
burst.  Figure  IV. 13a  shows  the  heterodyned  signal  from  WKU  and  MAIPU 
on  lines  one  and  two,  respectively.   Line  three  is  the  cross  correla- 
tion function  plotted  versus  time  for  a  shift  of  38  data  points,  and 
line  four  is  the  normalization.  Figure  IV. 13b  shows  the  cross 
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correlation  function  plotted  versus  time  for  other  shifts.  The 
shift  is  37,  39,  and  40  data  points  in  lines  one,  two  and  three, 
respectively.  The  integration  time  for  calculating  C  was  2.1  milli- 
seconds. It  should  be  obvious  from  the  two  parts  of  Figure  IV. 13 
that  the  cross  correlation  function  is  very  sensitive  to  relative 
shift.  The  value  of  best  shift  of  38  data  points  indicates  that  the 
clock  at  MAIPU  must  be  advanced  by  2.65  milliseconds  from  the  time 
that  was  indicated  by  the  corrections  given  to  BADE.  Figure  IV. 14 
is  an  illustration  of  the  full-wave  detected  data  from  WKU  and  MAIPU 
arranged  in  the  same  way  as  for  Figure  IV. 13.   Figure  IV. 14b  shows 
the  dependence  of  the  cross  correlation  function  on  relative  shift. 
The  deviation  from  the  value  for  the  best  shift  is  not  as  rapid  as 
for  the  heterodyned  data  and  is  the  chief  error  introduced  in  the 
technique.  Table  IV. 4  is  a  summary  of  the  best  relative  shifts  for 
the  data  from  the  three  pairs  of  stations. 


Stations 

WKU  -  FPC 
FPC  -  MAIPU 
WKU  -  MAIPU 


Table  IV. 4 

Best  Shift 
Heterodyned  Data     Detected  Data 
1.05  msec  1.05  msec 

1.54  msec         1.68  msec 
2.66  msec         2.66  msec 


Data  Pt. 

1200 

1200 

840 


Figure  IV. 15a  is  a  graph  of  the  dependence  of  the  cross  correla- 
tion function  on  the  relative  shift  for  the  entire  burst  as  received 
at  WKU  and  MAIPU.  The  integration  time  for  calculating  C2  was  112 
milliseconds,  which  covers  the  entire  burst  in  one  summation.  While  the 
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Figure  IV. 15 


—  Cross  Correlation  Measures  Plotted  as  Function  of 
Relative  Shift  for  the  WKU-MAIPU  Baseline,   (a)  C2 
for  Detected  Data  with  N=112  Milliseconds,   (b)  and 
(d)  C2  for  Detected  Data  in  the  Vicinity  of  Data  Points 
840  and  1240,  Respectively,  with  N=15  Milliseconds,  (c) 
and  (e)  Amplitude  of  the  Cross  Correlation  Function 
Near  the  Same  Two  Locations  with  N=2.1  Milliseconds. 


79 


curve  is  not  sharply  peaked,  it  is  obvious  that  the  best  shift  for  the 
entire  burst  is  within  400  microseconds  of  the  best  shift  for  the 
separate  sections  of  the  burst.  Figures  IV. 15b  and  c  are  graphs  of 
C  and  |c|  for  the  detected  and  heterodyned  data,  respectively,  near 
the  data  point  850.  Parts  d  and  e  are  for  C2  and  |c|  near  the  data 
point  1240  in  the  bursts  received  at  WKU  and  MAIPU.  The  time  axis  is 
calibrated  in  multiples  of  data  points  (1  dp  is  0.07  milliseconds) 
relative  to  the  best  shift  for  the  heterodyned  data  (2.66  milliseconds) 
All  curves  show  a  peak  in  the  vicinity  of  a  relative  best  shift  of 
zero.  Since  there  is  no  consistent  deviation  between  the  best  shift 
for  the  heterodyned  data  and  that  for  the  detected  data,  it  is 
concluded  that  there  is  no  measurable  beam  sweeping  effect,  and  that 
the  observed  slight  deviation  is  experimental  scatter  due  to  contamina- 
tion of  the  Jupiter  signal  by  uncorrelated  galactic  noise.  The 
uncertainty  in  timing  is  on  the  order  of  400  microseconds.   If  a  beam 
had  swept  past  the  two  stations  at  a  rate  slower  than  this,  its 
presence  would  have  been  detected. 

On  2  January  1969  the  distance  from  the  earth  to  Jupiter  was 
7.81  x  10  km.  The  north-south  baseline  in  equatorial  coordinates  for 
WKU  -  MAIPU  was  1.9  seconds  of  arc  as  seen  from  Jupiter.  The  east-west 
baseline  in  the  same  coordinate  system  was  0.50  seconds  of  arc.  The 
maximum  detectable  beam  sweep  rate  in  the  north-south  direction  was 
2°/second  and  was  0?5/second  in  the  east-west  direction.  For  a  beam 
sweeping  parallel  to  the  baseline,  the  maximum  detectable  sweep  rate 
was  on  the  order  of  2?l/second. 

The  theory  of  Goldreich  and  Lynden-Bell  contains  four  possible 
mechanisms  that  might  cause  a  beam  of  radiation  to  be  swept  across 
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the  earth.  Jupiter's  ionosphere,  where  the  lower  terminus  of  a  flux 
tube  is  located,  rotates  with  the  planet.  If  the  source,  moving  with 
the  ionosphere,  swept  the  beam  of  radiation  with  it,  the  beam  would 
rotate  with  an  angular  velocity  of  0?01  per  second.  The  flux  tube 
passes  through  Io  and,  in  this  theory,  moves  with  the  satellite. 
The  orbital  angular  velocity  of  Io  is  0?0023  per  second,  which  would 
be  the  sweep  rate  of  a  beam  influenced  by  this  satellite.  The  theory 
proposes  that  as  ionization  conditions  in  the  ionosphere  change  from 
those  at  night  to  the  daytime  state,  that  the  foot  of  the  flux  tube 
advances  from  a  position  below  Io  to  one  approximately  15°  ahead.  If 
the  beam  is  swept  by  this  mechanism,  its  angular  velocity  would  be  on 
the  order  of  0?004  per  second.  The  advance  is  assumed  to  occur  in 
about  20  minutes.   It  is  obvious  that  any  of  these  beam  sweep  rates 
would  be  easily  detectable  by  the  method  proposed  above.  A  group  of 
electrons  spirals  along  a  field  line  that  lies  between  Io  and  the 
ionosphere  of  Jupiter.  If  it  is  assumed  that  the  radius  of  curvature 
of  the  field  line  in  this  region  is  on  the  order  of  5  times  the  radius 
of  the  planet  and  that  the  electrons  are  moving  with  a  velocity  of 
0.1  times  the  speed  of  light,  a  beam  of  radiation  generated  by  them 
would  sweep  at  a  rate  on  the  order  of  4?8  per  second.  This  sweep  rate 
is  just  beyond  the  range  of  measurement  by  the  data  that  has  been 
analyzed  in  this  chapter.  The  three  sweep  mechanisms  mentioned  earlier 
are  definitely  ruled  out,  however. 

Douglas  and  Smith  (1967)  proposed  a  theory  that  accounts  for  the 
structure  of  the  L-Burst.  They  attribute  the  L-Burst  to  an  isophotal 
pattern,  caused  by  inhomogeneities  in  the  solar  wind,  drifting  across 
the  receiving  stations.  They  support  their  theory  with  convincing 
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evidence  based  upon  the  time  of  arrival  of  a  given  L-Burst  at  several 
receiving  stations.  Their  model  predicts  that  the  easternmost 
station  on  a  baseline  parallel  to  the  ecliptic  will  receive  the  burst 
before  any  stations  to  the  west  if  Jupiter  is  before  opposition. 
The  westernmost  station  will  receive  the  burst  earliest  if  Jupiter  is 
after  opposition.  The  L-Burst  analyzed  in  this  chapter  was  received 
when  Jupiter  was  82  days  before  opposition.  Their  data  shows  that 
for  a  baseline  length  of  about  100  km  that  the  signal  will  arrive  at 
the  station  to  the  east  about  0.1  second  before  it  arrives  at  the 
western  station  for  this  time  (82  days  before  opposition) .  The  WKU  - 
MAIPU  baseline,  when  projected  on  the  ecliptic,  B    has  a  length  of 
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4580  km  on  a  plane  perpendicular  to  the  line-of-sight  from  the  earth 
to  Jupiter.  By  scaling  the  time  delay  in  a  linear  fashion  proportional 
to  the  relative  baseline  lengths,  it  is  seen  that  the  expected  delay 
time  should  be  on  the  order  of  4.5  seconds.  No  such  effect  was 
observed  in  this  case.   It  is  therefore  concluded  that  the  source  itself 
may  exhibit  intensity  variations  that  are  at  least  as  long  as  the 
burst  that  has  been  analyzed  here.  Although  most  L-Burst  envelopes  are 
probably  shaped  by  the  solar  wind,  some  must  be  intrinsic  to  the  source 
itself. 

IV-5.  The  Stability  of  the  Position  of  the  Source 

While  the  phase  interferometer  can  give  no  information  about  the 
size  of  the  source  that  radiates  coherently,  the  fringes  are  a  good 
source  of  information  about  the  relative  motion  of  the  source  during 
the  time  of  emission.  Examination  of  Figure  IV. 4  will  reveal  no 
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consistent  departure  of  the  fringe  rate  or  rapid  changes  of  phase 
from  that  shown  at  the  beginning  of  the  burst.  There  are  numerous 
apparent  losses  of  coherence,  but  these  can  be  correlated  with 
intervals  of  low  signal  level.  The  effect  of  these  low  levels  is  to 
allow  the  uncorrelated  noise  of  the  galactic  background  to  cause 
temporary  loss  of  coherence.  When  the  signal  increases  again,  the 
fringe  regains  its  coherence,  with  no  shift  in  phase.  The  inter- 
pretation of  this  effect  is  that  the  source  remained  in  essentially  the 
same  position  in  the  interferometer  antenna  pattern  during  the 
emission  of  the  radiation.  The  spatial  resolution  is  on  the  order 
of  0.1  part  of  a  fringe  or  about  0.05  seconds  of  arc.  This  corresponds 
to  the  source  remaining  within  245  km  of  its  original  position  during 
the  burst. 


CHAPTER  V 
ANALYSIS  OF  DATA  RECEIVED  IN  1970 

V-l.  Instrumentation,  Timing  and  Baseline  Orientation 

In  1970  the  receivers  that  are  diagrammed  in  Figures  III. 2,  III. 3, 
and  III. 4  were  put  on  the  air  at  the  University  of  Florida  Radio 
Observatory,  Old  Town,  Florida,  and  at  the  Observatorio  Radioastronomico 
de  Maipu  of  the  Universidad  de  Chile  near  Maipu,  Chile.  The  geographical 
coordinates  of  these  stations  were  given  in  Chapter  IV.  The  great 
circle  distance  between  UFRO  and  MAIPU  is  7384  km  and  the  chordal 
distance,  D,  is  6980  km  (419,000  A);  the  angle  between  first  nulls  for 
the  interferometer  is  0.49  seconds  of  arc.   Figure  V.la  shows  a  view 
of  the  earth  as  seen  from  the  North  Pole  of  the  ecliptic  on  30  April 
1970  at  the  time  the  bursts  that  are  analyzed  in  this  chapter  were 
received.  The  celestial  coordinates  of  Jupiter  on  that  date  were  as 
follows:  right  ascension  =  13  53m54s  and  declination  =  -10o5'41". 
The  right  ascension  of  the  sun  was  2  27m22s. 

The  various  baseline  projections  that  will  be  used  in  the  interpreta- 
tion of  the  data  are  shown  in  Figure  V.lb.  The  projections  are  drawn 
in  a  plane  that  is  perpendicular  to  the  line-of-sight  from  the  stations 
to  Jupiter.  D   is  the  component  of  D  parallel  to  the  polar  axis  of 
the  earth  projected  on  the  plane  perpendicular  to  the  line-of-sight 
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Figure  V.l  --  (a)  View  of  Earth,  Jupiter  and  the  Sun  from  Above 

the  North  Pole  on  30  April  1970,  (b)  D   and  D 

NS      gw, 
Projected  on  the  Plane  Perpendicular  to  the  Line- 
of-Sight  to  Jupiter. 
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to  Jupiter.  D   is  the  component  of  D  parallel  to  the  earth's  equatorial 
plane  projected  on  the  same  plane.   Angle  a  lies  between  D   and  the 

No 

North  Pole  of  the  ecliptic  in  the  same  plane  and  had  a  value  of  20?5 

when  the  data  was  received.  The  various  baseline  components  had  the 

following  lengths:  D   =  6560  km  and  Dcir  =  205  km,  with  MAIPU  east  of 

No  cW 

UFRO  using  the  directions  for  the  side  of  the  earth  toward  Jupiter. 

To  insure  the  accurate  calibration  of  the  locally  generated  time 
signal  at  MAIPU  with  respect  to  real  time,  a  rubidium  vapor  clock  was 
provided  on  loan  from  the  National  Radio  Astronomy  Observatory.  Without 
interrupting  its  operation,  the  clock  could  be  transported  to  the  NASA 
satellite  tracking  station  near  Santiago  for  calibration.  This  procedure 
was  performed  at  biweekly  intervals  during  the  observing  season.  The 
data  channel  was  multiplexed  in  such  a  way  that  the  WWV  voice  announce- 
ment was  recorded  for  20  seconds  at  5-minute  intervals  and  the  rubidium 
vapor  clock  signal  was  recorded  for  three  seconds  every  minute  on 
the  minute. 

While  it  was  hoped  to  use  signals  transmitted  by  the  East  Coast 
LORAN-C  chain  to  provide  very  accurate  timing  at  UFRO,  the  audio 
preamplifier  in  the  LORAN-C  receiver  was  found  to  be  saturating,  and 
the  station  identifiers  could  not  be  decoded  on  chart  recordings  made 
from  slowed-down  magnetic  tapes.  Local  timing  was  determined  as 
accurately  as  possible  by  correcting  the  WWV  time  marks  for  transmission 
time  from  Boulder,  Colorado,  which  amounted  to  8.3  milliseconds  (Hewlett- 
Packard  Company,  1965) . 

The  data  format,  as  recorded  on  magnetic  tape  at  UFRO  and  MAIPU, 
was  given  in  Chapter  III.  Just  prior  to  the  installation  of  the 
receivers  at  both  stations,  a  fourth  timing  signal  was  added  to  the 
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local  time  channel  for  the  purpose  of  synchronizing  an  analog-to- 
digital  converter.  The  sync  pulse  had  a  frequency  of  15  kHz  and  was 
derived  by  a  separate  divide  chain  from  the  1.8  MHz  crystal  at  each 
station.  Synchronization  of  the  A-to-D  converter  to  this  signal  caused 
the  resulting  digitization  rate  to  be  constant  with  respect  to  the 
local  time  standard  (hence  constant  in  real  time) .  Since  this 
effectively  eliminated  the  variation  in  digitization  rate  caused  by  tape 
stretch  and  speed  instabilities  in  the  magnetic  tape  drives,  the  most 
complex  part  of  the  program,  BADE,  could  be  eliminated. 

V-2.  The  Nature  of  the  Received  Data 

Figures  V.2  and  V.3  reproduce  chart  recordings  of  the  envelope 
detected  data  as  received  at  UFRO  (both  channels)  and  MAIPU  beginning 
at  9h35m51?5  U.T.  on  30  April  1970.  The  time  covered  by  the  chart 
records  is  about  4.5  seconds.  The  taller  marks  in  the  time  channels, 
which  are  shown  on  alternate  lines  just  below  the  corresponding  data 
channels,  occur  at  one-second  intervals,  and  the  shorter  marks  occur 
at  20-millisecond  intervals.  The  one-millisecond  marks  and  the  15  kHz 
sync  tone  are  not  distinguishable. 

There  were  approximately  68  S-Bursts  in  the  time  period  beginning 
at  9  35  51?6  U.T.  and  ending  six  seconds  later.  Of  these  68  bursts, 
only  7%  were  present  at  one  station  and  missing  at  the  other.  This 
analysis  will  be  restricted  to  the  eleven  bursts  that  are  indicated 
in  Figures  V.2  and  V.3.  Table  V.l  gives  the  beginning  times  of  the 
intervals  that  were  digitized  in  terms  of  the  locally  generated  times 
at  each  station  and  in  terms  of  U.T.  The  indicated  times  are  in 
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seconds  after  9h35m51s.  The  digitized  period  extended  for  1920  data 
points  C51.3  milliseconds)  from  the  times  given.  All  bursts  that  are 
analyzed  are  contained  within  these  intervals. 


Table  V.l 


Burst  No. 

MAIPU   (local 

time) 

UFRO   (local 

time) 

U.T. 

1 

0.618 

0.610 

0.6193 

2 

1.168 

1.160 

1.1693 

3 

1.288 

1.280 

1.2893 

4 

2.998 

2.990 

2.9993 

5 

3.208 

3.200 

3.2093 

6 

3.778 

3.770 

3.7793 

7 

3.928 

3.920 

3.9293 

8 

4.078 

4.070 

4.0793 

9 

4.138 

4.130 

4.1393 

10 

4.213 

4.205 

4.2143 

11 

4.7607 

4.750 

4.7619 

Two  types  of  bursts  were  selected  for  analysis.  The  first  type, 
consisting  of  numbers  1,  3,  4  and  9,  were  of  short  duration  (about 
20  milliseconds)  and  had  sharp  rise-times  and  great  intensity.  The 
second  type,  numbers  2,  5,  6,  7,  8,  10  and  11,  were  of  longer  duration 
(about  50  milliseconds)  but  are  still  classified  as  S-Bursts  according 
to  the  system  of  Gallet  mentioned  in  Chapter  I.  The  ratio  of  Jupiter 
signal  +  galactic  noise  to  galactic  noise  alone  was  exceptionally  high 
for  all  of  these  bursts  as  compared  with  that  for  the  L-Burst  discussed 
in  Chapter  IV. 
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Figure  V.4  reproduces  chart  recordings  and  incremental  plotter 
graphs  of  burst  number  5.  The  chart  recordings  are  shown  for  MAIPU 
data  and  time  channels  on  line  (b)  and  for  UFRO  data  and  time  channels 
on  line  (d) .  The  taller  marks  visible  in  the  time  channels  are  20- 
millisecond  marks,  and  the  shorter  are  1-millisecond  marks.  The 
computer  incremental  plotter  graphs  are  shown  for  MAIPU  on  line  (a) 
and  for  UFRO  on  line  (c) .  The  time  scale  for  the  computer  plots  is 
expanded  by  about  eight-to-one  over  that  for  the  chart  recordings. 
The  digitization  rate  used  for  this  data  was  37,500  data  points/real- 
time second.  The  feature  marked  by  A  on  all  four  lines  corresponds  to 
the  same  feature  in  each  burst. 

V-3.  Analysis  of  the  Digitized  Data 

IBM  cards  containing  1920  data  points  for  each  burst  from  each 
station  served  to  transport  the  data  from  the  IBM  1800  at  the  Univer- 
sity of  Florida  Medical  Center,  where  it  was  digitized,  to  the  IBM 
360/65  at  the  University  of  Florida  Computing  Center,  where  it  was 
processed.  The  two  phase-quadrature  channels,  which  will  be  referred 
to  as  UFRO-A  and  UFRO-B,  were  digitized  separately.  The  Appendix 
contains  a  listing  of  the  changes  made  to  CROCO  in  order  to  process 
this  data.  The  principal  changes  were  necessary  because  of  the 
different  number  of  data  points  in  a  burst  as  compared  with  the 
L-Burst  in  Chapter  IV.  A  program  was  written  in  IBM  360  Assembler 
Language  to  convert  the  column  binary  formatted  data  cards  to  files 
on  magnetic  tape  for  simpler  access  by  CROCO.  The  data  for  a  given 
burst  from  a  station  was  arranged  on  the  tape  in  the  form  of  24  records 
of  80  data  points  each.  The  first  station  written  for  each  burst 
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was  MAIPU,  the  second  was  UFRO-A  and  the  third  was  UFRO-B.  The  first 
several  statements  in  the  adapted  version  of  CROCO  select  the  appropriate 
burst  using  IFIL1  (which  ranges  from  1  through  11)  and  read  the  data 
for  that  burst  into  arrays  XI  (which  receives  MAIPU  data)  and  IGV 
(which  receives  both  channels  of  UFRO  data) .  IFIL2  selects  which 
channel  of  UFRO  data  will  be  cross  correlated  with  MAIPU  data  by  CROCO. 
The  baseline  adjustment  portion  of  BADE  was  added  to  CROCO  in  the 
vicinity  of  statements  104  and  501  since  the  phase- locking  of  the 
A-to-D  conversion  procedure  had  eliminated  the  need  for  the  remainder 
of  the  program  BADE. 

In  certain  cases  it  was  necessary  to  provide  detected  data  to 
CROCO.  Since  only  heterodyned  data  was  digitized,  several  statements 
were  placed  in  the  input  portion  of  CROCO  to  square- law  detect  the 
heterodyned  data,  filter  it  using  a  1-millisecond  time  constant, 
and  adjust  its  baseline  to  the  average  of  the  burst.  The  added 
sections,  which  are  shown  in  the  Appendix,  could  easily  be  removed 
if  it  were  desired  to  return  CROCO  to  the  form  used  to  process 
heterodyned  data. 

Two  types  of  analyses  were  performed  on  the  S-Burst  data.  The 
bursts  that  were  of  short  duration  and  extreme  intensity  were  used  to 
test  the  sweeping  beam  hypothesis.  The  longer  bursts  were  used  in  an 
effort  to  determine  the  fringe  rate  of  the  interferometer.  In  order 
to  perform  the  sweeping  beam  test  it  was  necessary  to  cross 
correlate  both  the  heterodyned  and  square- law  detected  data.  The 
fringe  rate  and  fringe  stability  calculations  could  be  completed  using 
only  the  cross  correlated  heterodyned  data. 
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V-4.  Burst  Arrival  Time  Analysis 

A  method  to  improve  the  timing  precision  using  the  properties  of 
the  Jovian  radio  burst  was  demonstrated  in  Chapter  IV.  The  analysis, 
using  1970  data  with  its  better  timing  and  faster  rise- times,  was 
expected  to  raise  the  maximum  detectable  sweep  rate  over  that  found 
using  the  burst  in  Chapter  IV  due  to  the  much  greater  timing  resolution 
used  in  digitizing  the  later  data  (37,500  data  points/second  as 
compared  with  14,200  data  points/second  for  the  1969  data) . 

Figures  V.5,  V.6,  V.7,  and  V.8  show  the  various  cross  correlation 

functions  plotted  against  relative  shift,  t,  for  the  short  bursts 

(1,  3,  4  and  9).  Since  these  bursts  have  such  fast  rise-times,  it  was 

felt  that  the  best  shifts  determined  for  the  envelopes  would  be  more 

reliable  than  those  found  in  Chapter  IV  for  the  1969  L-Burst.  The  (a) 

line  of  each  of  the  four  figures  represents  a  plot  of  C2,  the  output 

of  CROCO  using  detected  da, a,  versus  the  relative  shift,  t,  for  an 

integration  interval  of  18.7  milliseconds  (700  data  points),  which 

spans  each  burst  completely.  The  relative  shift  means  that  data  point 

E    in  Equation  4)  (Chapter  IIT)  will  be  multiplied  by  E'    ,  where 
Oj  0  j+k 

k  =  (digitization  rate  )X(  t) .  Here,  station  1  (as.  indicated  by  the 
subscript  onEQ)was  MAIPU  and  station  2  was  either  UFRO-A  or  UFRO-B. 
In  the  figures  it  is  seen  that  the  range  of  shifts  indicates  that 
the  MAIPU  data  table  is  being  read  10  to  40  data  points  ahead  of  the 
UFRO  table.  The  (b)  lines  of  Figures  V.5,  V.6,  V.7  and  V.8  show 


/C^  +  C2  plotted  versus  relative  shift,  where  C  is  the  value  of  C 

for  a  time  near  the  center  of  the  burst  using  heterodyned  data  from 

MAIPU  and  UFRO-A  and  C  is  similar  except  that  data  from  UFRO-B  is 
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used  instead  of  UFRO-A.  The  integration  interval,  N,  is  6.6  milli- 
seconds (250  data  points)  and  is  located  near  the  most  intense  portion 
of  the  burst. 

It  was  shown  in  Chapter  II  that  the  best  shift,  i.e.,  that  value 
of  t  that  maximized  this  square-root,  could  be  found  regardless  of  the 
phase  of  the  fringes.  In  a  burst  as  short  as  the  ones  under  considera- 
tion, it  is  very  unlikely  that  the  fringe  would  be  passing  through  a 
maximum  during  the  burst,  making  this  technique  necessary.  The 
oscillatory  pattern  of  maxima  and  minima  is  due  to  the  small  ratio 
of  bandwidth  to  center  frequency  that  was  used.  The  receivers  were 
operating  with  Collins  mechanical  filters  having  2.1  kHz  bandwidths  with 
passbands  centered  at  5  kHz.  The  Jupiter  signal,  as  shown  in  the 
computer  plots  in  Figure  V.4a  and  c,  is  seen  to  be  quasi-sinusoidal 
with  the  frequency  varying  in  a  random  manner  but  over  a  relatively  small 
range  as  restricted  by  the  mechanical  filters  in  the  receiver  I.F. 
amplifiers.   It  may  be  recalled  that  the  autocorrelation  function  for 
a  sine  wave  consists  of  a  similar  oscillation,  except  that  all  of  the 
maxima  are  the  same  height  and  all  of  the  minima  are  zero.  A  given 
maximum  and  the  adjacent  minimum  are  one-half  of  the  period  of  the  sine 
wave  apart.  In  Figure  V.5b  it  is  seen  that  the  maxima  vary  in  height 
in  an  orderly  way  such  that  a  central  maximum  can  be  detected  at  a 
relative  shift  of  25.5  data  points.  The  spacing  between  maxima 
averages  about  3.7  data  points  (1  millisecond)  which  is  one-^half  of 
the  period  of  the  center  frequency  of  the  passband  (5  kHz) .  This 
apparent  halving  of  the  "period"  of  C  versus  t  as  compared  to  the 
autocorrelation  example  above  is  caused  by  the  plotting  of 


^CA  +  CB  which  serves  to  "rectify"  the  cross  correlation  function. 
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Table  V.2  summarizes,  the  values  of  the  best  shift,  t    ,  for 

best' 

both  heterodyned  and  detected  data  for  the  bursts  shown  in  Figures 
V.5,  V.6,  V.7  and  V.8. 


Table  V.2 

Burst 

No. 

T 

best 
(het.  data) 
Data  Points 

25.5 

x 
best 

(det.  data) 

Data  Points 

28.0 

At  Data 
Points 

-2.5 

Max.  Timing 
Error 

1 

(±1.8  d.p.) 

3 

29.4 

27.6 

+  1.8 

(+1.8  d.p.) 

4 

27.2 

25.6 

+  1.6 

(±1.8  d.p.) 

9 

25.4 

25.0 

+0.4 

(±1.8  d.p.) 

The  quantity,  At,  is  calculated  from  At  =  t    (het.  data)  - 

best 

TbestCdet'  data)- 

There  are  two  mechanisms  that  must  be  considered  in  order  to 

determine  the  cause  of  the  observed  delays,  At.  The  first  concerns  the 

accuracy  in  timing  the  arrival  of  a  pulse  considering  that  the  pulse 

transversed  a  filter  in  a  receiver  which  had  a  bandwidth  of  2  kHz.  The 

time  for  a  significant  change  in  signal  level  allowed  by  such  a  filter 

is  on  the  order  of  500  microseconds.  However,  the  envelope  cross 

correlation  technique  should  be  able  to  detect  the  arrival  times  at  least 

an  order  of  magnitude  more  accurately  than  this  since  the  technique 

consists  of  matching  the  slopes  of  the  rise-  and  fall-times  of  bursts 

received  at  different  places.  Any  delay,  At,  less  than  50  microseconds 

(1.8  data  points)  should  be  interpreted  as  being  on  the  edge  of,  or 

beyond,  the  limit  of  detectability  for  data  received  using  the  2.1  kHz 

bandwidth.  The  delay,  At,  measured  for  burst  number  9  is  clearly  below 
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the  minimum  detectable  delay  since  it  has  a  value  of  0.4  data  points 
(11  microseconds).  On  this  account,  the  evidence  from  burst  9  will 
not  be  considered  further. 

The  second  mechanism  for  consideration  involves  a  beam  of  radiation 
from  Jupiter  being  swept  across  the  earth  so  that  it  is  received  at  one 
station  a  short  time  before  it  is  received  at  another  station  some 
distance  away.  Using  the  example  given  in  Figure  IV. 12,  the  delay, 
At,  would  have  a  value  of  -3  time  units  (not  data  points)  if  station  1 
were  interpreted  as  UFRO  and  station  2  as  MAIPU.  Notice  that  the  best 
relative  shift  for  the  heterodyned  data  in  the  example  is  zero  time 
units  and  is  +3  time  units  for  the  detected  data  (the  envelopes),  i.e., 
the  MAIPU  data  table  is  read  about  3  time  units  later  than  the  UFRO 
data  table  in  calculating  the  best  shifted  cross  correlation  function 
for  the  detected  data.  Therefore,  using  the  example  as  a  model,  if  the 
delays  were  caused  by  a  sweeping  beam,  those  At's  with  a  negative  sign 
would  indicate  that  the  beam  had  a  component  sweeping  from  UFRO  to 
MAIPU  and  a  positive  sign  would  indicate  a  sweep  in  the  reverse  direc- 
tion. Table  V.3  shows  the  values  of  At  expressed  in  milliseconds  for 
bursts  1,  3  and  4  and  the  angular  sweep  rate  parallel  to  the  baseline 
components,  D,  DNS  and  D  ,  These  components  subtended  angles  of 
2.18  seconds  of  arc,  1.05  seconds  of  arc  and  0.064  seconds  of  arc, 
respectively,  as  seen  from  Jupiter  at  the  time  the  bursts  were  received. 
The  At's  for  bursts  3  and  4  are  at  the  time  resolution  limit  of  the 
system,  while  that  for  burst  1  is  just  above  the  limit  (within  the  range 
of  detection).  They  are  shown  to  indicate  that,  if  sweeping  beams  are 
present,  their  sweep  rates  must  have  been  at  this  rate  or  greater. 
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Tabl 

e  V.3 

Burst  No. 

At (msec) 

Parallel 

Sweep 
to  D 

Minimum  Possible 
Rate  Cdegrees/second) 

Parallel  to  DXT„   Para] 

NS 

8.7 

•  lei  to  DEW 

1 

-0.065 

9.3 

0.27 

3 

+0.048 

12.5 

11.8 

0.37 

4 

+0.043 

14.0 

13.1 

0.41 

Max.  Det. 

±0.050 

11.9 

11.1 

0.35 

In  Chapter  IV  it  was  indicated  that  the  maximum  east-west  sweep  rate 
parallel  to  the  ecliptic,  as  caused  by  Jupiter's  rotation,  was  on  the 
order  of  0?01/second.  It  is  obvious  from  Table  V.3  that  no  beams  were 
detected  with  such  a  slow  sweep  rate.  A  possible  north-south  sweeping 
might  have  been  contributed  by  the  motion  of  bunches  of  electrons  along 
the  curved  magnetic  field  lines  as  they  radiated.  The  model  assumed  in 
Chapter  IV  predicted  sweep  rates  on  the  order  of  5°/second  or  higher. 
It  can  be  concluded  from  the  data  presented  in  Table  V.3  that  the  bursts 
analyzed  in  this  section  were  not  connected  with  beams  that  swept  at  a 
rate  less  than  9°/second,  if  the  beams  were  swept  at  all.  Since  the 
At  for  burst  1  was  within  the  time  resolution  limit,  it  should  be 
interpreted  as  a  reason  to  make  more  observations  with  broader  bandwidth 
equipment.  This  one  data  point  does  not  establish  the  existence  of 
sweeping  beams,  but  it  should  encourage  an  extended  search. 
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V-5.  Measurement  of  the  Size  of  the  Source 

Figures  V.9  and  V.10  complete  the  graphical  presentation  of  the 
cross  correlation  function  plotted  against  relative  shift  for  the  rest 
of  the  bursts.  Table  V.4  is  a  summary  of  the  best  relative  shifts 
for  these  bursts. 

Table  V.4 
Burst  No.  Best  Relative  Shift  (data  points) 

2  29.6 

5  27.0 

6  23.5 

7  28.2 

10  23.5 

11  24.3 

Burst  number  8  was  not  included  since  an  error  in  locating  the 
digitization  interval  for  UFRO-B  prevented  using  that  data  to 
calculate  C  . 

There  are  two  features  in  the  figures  that  must  be  considered  at 
this  point.  The  first  is  the  failure  of  each  burst  to  have  the  same 
best  relative  shift.  This  accounted  for  by  considering  the  physical 
relationship  of  the  A-to-D  converter  sync  tone  and  the  data  channel  on 
the  analog  magnetic  tape.  A  simple  calculation  shows  that  a  dynamic 
skew  of  0?6  (a  line  drawn  perpendicular  to  the  edge  of  the  tape 
changes  by  this  angle  with  respect  to  the  head  gap  as  the  tape  moves) 
can  account  for  the  6  data  point  range  in  best  relative  shift.  This 
calculation  assumes  that  the  track  separation  is  0.125  inches  and  the 
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A3 


A       B 


/C2*  C2 
A       B 


Figure  V.9  —  /C*  ♦  C|     Versus  Relative  Shift  for  (a)  Burst  Number  2, 
(b)   Burst  Number  5  and   (c)   Burst  Number  6. 
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20  22         24       26         28  30       32 

RELATIVE  SHIFT  (data  points) 
/C2+c2  versus  Relative  Shift  for  (a)  Burst  Nunber  7, 

A        D 

(b)  Burst  Number  10  and  (c)  Burst  Number  11. 
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recording  rate  is  7.5  inches/second.  This  value  of  dynamic  skew  has 
been  detected  in  the  0.25  inch  magnetic  tape  recorders  used  for  recording 
the  data.  In  the  future  the  sync  tone  will  be  placed  in  the  data  channel 
to  eliminate  this  problem. 

The  other  feature  that  is  obvious  in  Figure  V.5,  for  example,  is 


that  the  minima  in  /C2  +  C2  versus  relative  shift  do  not  go  to  zero. 

A    B 

If  the  local  oscillators  at  UFRO  that  provided  the  L.O.  injection  for 

the  A  and  B  channels  had  been  exactly  90°  out  of  phase,  these  minima 

would  have  been  zero  and  the  peak  maximum  would  have  been  unity.  The 

phase  angle  can  be  calculated  for  a  burst  that  has  a  negligible  amount 

of  galactic  contamination  and  can  be  compared  with  other  similarly 

strong  bursts  to  corroborate  the  possibility  of  this  error.  Using  the 

equations  on  page  31  in  Chapter  II  and  choosing  suitable  values  for 

the  variables,  the  error  can  be  evaluated  as  shown  in  the  next  paragraph. 

Let  v,       -   C,  =  Vn(s.)  and  V„  ,   .   =  C  =  V„ (s  )  cos  (90°  -  a), 
1st  pair    A    0     X  2nd  pair    B    0 v  A      v      J ' 

where  a   is  the  deviation  in  the  relative  phase  at  UFRO  from  90°.  Solving 
the  equations  simultaneously  for  a   gives  the  following: 


sin  o 


c2  +  c2 

A          B 

-    1 

_VS  W 

Using  burst  1,  which  is  relatively  strong 


and  is  therefore  more  likely  to  have  a  fringe  visibility  close  to 
unity,  o  is  found  to  be  34°.  A  comparison  of  the  value  of 


^CA  +  CB  for  t^ie  best  relative  shift  from  other  strong  bursts  shows 
that  0  does  not  vary  appreciably  over  the  4  seconds  covered  by  the 
data  analysis.  This  is  as  would  be  expected  for  an  equipment  problem 
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such  as  this .  It  is  not  possible  with  the  present  data  to  determine 
the  value  of  the  fringe  visibility  without  assuming  that  it  is  unity 
for  at  least  one  burst.  Since  CA  does  reach  a  value  of  0.8  for  burst  1, 
this  is  not  a  bad  assumption  considering  that  C.  plotted  as  a  function 
of  time  may  not  be  at  a  fringe  maximum  during  burst  1.  Using  this 

assumption,  the  o-corrected  value  of  /cf  +  C*  averaged  over  all  bursts 

n  B 

(taking  into  account  the  varying  degree  of  galactic  noise  contamination) 
is  1.01  ±  0.07.  The  maximum  angular  diameter  of  an  incoherent, 
Gaussian  distributed  source  is  0.1  seconds  of  arc  C320  km  at  Jupiter). 
The  mean  angular  size  is  0.0  seconds  of  arc,  with  an  error  of  +0.1 
seconds  and  -0.0  seconds. 

It  was  hoped  at  the  time  that  the  data  was  digitized  that  the 
fringe  rate  (CA  plotted  versus  time)  would  be  such  that  data  spaced 
over  the  approximately  4-second  interval  chosen  would  allow  the 
reconstruction  of  the  fringes.  A  graph  of  C  versus  time  is  given  in 

Figure  V.ll,  where  the  bursts  and  their  values  of  C  are  labelled 

A 

Some  segments  of  data  were  long  enough  to  allow  a  portion  of  the  fringe 

to  be  seen.  In  this  case  lines  were  drawn  to  show  the  part  of  the 

fringe  that  is  visible.  It  was,  in  general,  impossible  to  determine 

one  curve  that  fit  all  of  the  values  of  C  plotted  versus  time.  While 

A 
this  might  be  interpreted  as  due  to  the  "jumping  about"  of  the  source 

causing  the  fringes  to  fail  to  maintain  phase,  it  is  felt  that  due  to 

the  lack  of  values  of  C  between  those  plotted,  the  interpretation 

should  be  withheld  until  the  intermediate  data  points  can  be  calculated. 

Using  the  fraction  of  the  fringe  period  that  is  present  in  burst  7, 

a  value  for  the  fringe  rate  at  that  time  can  be  calculated.  It  is  on  the 
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order  of  9  fringes/second,  which  compares  favorably  with  the  predicted 
value  of  9.3  fringes/second  due  to  the  earth's  rotation. 

V-6.  Brief  Summary  of  the  1970  Data  Analysis 

While  an  upper  limit  of  11.9  degrees/second  has  been  placed  on  the 
beam  sweep  detection  capability  of  the  system  used  in  1970  by  the  narrow 
bandwidth  of  the  receivers,  it  was  shown  that  at  least  one  burst  seemed 
to  indicate  a  sweep  rate  a  little  less  than  this  value,  about  9  degrees/ 
second.   It  is  concluded  that  while  all  of  the  apparent  sweep  rates 
are  so  close  to  the  limit  of  detectability  that  none  may  be  significant, 
there  is  sufficient  evidence  for  beam  sweeping  to  justify  a  continua- 
tion of  the  search  with  wider  bandwidth  receivers.  If  the  observed 
apparent  sweep  rates  are  real,  they  would  indicate  beaming  from  electron 
bunches  moving  along  curved  field  lines.  Sweep  rates  due  to  rotation 
effects  would  be  much  lower  and  are  therefore  ruled  out. 

The  maximum  fringe  amplitude  actually  measured  was  0.8,  but  it 
was  felt  that  this  value  corresponds  to  a  point  not  at  a  fringe  maximum. 
For  the  baseline  used,  the  a-corrected  /C2  +  c2  averaged  over  all  bursts 
having  a  value  of  1.01  ±  0.07  indicates  that  the  maximum  size  of  the 
source  (as  given  by  the  extreme  lower  limit  of  0.94)  is  less  than  0.1 
seconds  of  arc,  if  the  source  is  assumed  to  be  an  incoherent,  Gaussian 
distributed  radiator. 

Since  so  few  segments  of  the  fringe  were  calculated,  no  evidence 
could  be  cited  concerning  the  stability  of  the  fringes  and  hence  of  the 
source.  It  is  indicated  that  in  the  future,  when  A-to-D  conversion 
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facilities  become  available  that  remove  the  limit  on  the  number  of  data 
points  that  can  be  handled,  this  data  should  be  reprocessed  for  the 
purpose  of  examining  the  source  position  stability  problem. 


CHAPTER  VI 
ANALYSIS  OF  DATA  RECEIVED  IN  1971 

VI-1.  Instrumentation,  Timing  and  Baseline  Orientation 

The  receiving  systems  described  in  the  previous  chapter  for  the 
1970  Jupiter  season  were  used  again  in  1971  with  a  few  modifications. 
Two  more  stations  were  added  to  the  University  of  Florida  interferometer 
network.  A  receiver  identical  to  the  one  used  in  Chile  was  put  on  the 
air  at  Western  Kentucky  University  Radio  Observatory  near  Bowling 
Green,  Kentucky.  A  commercial  receiver  employing  the  same  local 
oscillator  frequency  synthesis  technique  as  used  in  the  other  inter- 
ferometer receivers  was  operated  by  members  of  Southwest  Research 
Institute  at  a  radio  observatory  near  San  Antonio,  Texas.  The 
observing  sessions  during  the  1971  apparition  produced  several  reels  of 
magnetic  tape  containing  mostly  L-Burst  activity.  A  single  L-Burst, 
found  on  tapes  from  UFRO  and  MAIPU,  is  discussed  in  this  chapter.  There 
were  several  promising  bursts  on  the  tape  from  SWRI  Radio  Observatory 
that  were  seen  on  chart  recordings  of  the  detected  data  from  that 
station.  No  chart  recordings  have  been  made  from  the  tapes  from  WKU. 

The  removal  of  the  rubidium  vapor  clock  that  was  used  at  MAIPU  in 
1970  required  that  alternate  means  of  calibrating  the  local  clock  be 
used.  The  one-second  ticks  generated  by  a  stable  crystal  oscillator 
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at  the  National  Observatory  in  Chile  were  used  to  key  a  tone  which 
amplitude  modulated  a  radio  transmitter.  The  transmitted  signal  was 
received  simultaneously  at  MAIPU  and  at  the  NASA  satellite  tracking 
station  near  Santiago,  Chile.  Since  the  distances  to  the  receiving 
stations  are  short  and  are  about  equal,  the  time  correction  supplied 
by  the  NASA  station  could  be  applied  directly  to  the  local  one-second 
ticks  generated  at  MAIPU.  The  NASA  time  correction  was  not  available 
when  the  data  was  cross  correlated  so  the  correction  will  be  compared 
with  the  value  of  the  best  relative  shift  presented  later  in  this  chapter. 

At  the  UFRO  the  LORAN-C  receiver,  having  been  repaired,  was  used 
to  supply  the  LORAN-C  bursts  to  the  tape  recorder.  The  one-second 
tick  transmitted  by  the  East  Coart  chain  master  station  at  Cape  Fear, 
South  Carolina,  was  used  to  calibrate  the  locally  generated  time  marks 
at  the  UFRO.  Since  the  LORAN-C  bursts  were  time  multiplexed  onto  the 
data  channels,  the  time  correction  was  applied  after  the  slowed-down 
tapes  were  displayed  on  chart  recordings.  Transmission  time  from  the 
master  station  was  taken  into  account.  The  method  used  is  outlined  by 
L.  D.  Shapiro  in  two  papers  (1968a  and  1968b). 

The  L-Burst  discussed  in  this  chapter  was  received  at  UFRO  and 
MAIPU  on  12  April  1971  at  8h51m23f387  U.T.  Figure  VI. la  shows  the 
directions  to  the  sun  and  to  Jupiter  as  seen  from  a  point  above  the 
earth  at  the  North  Pole  of  the  ecliptic.  The  figure  is  a  projection 
of  these  directions  in  the  plane  of  the  ecliptic.  The  right  ascension 
of  Jupiter  was  16h16m38s,  and  its  declination  was -20° 17 '30".  The 
right  ascension  of  the  sun  at  the  time  the  burst  was  received  was 
lh20m20s . 
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Figure  VI. 1  —  (a)  View  of  Earth,  Jupiter  and  the  Sun  from  Above 
the  North  Pole  on  12  April  1971,  (b)  DNS  and  D 

Projected  on  the  Plane  Perpendicular  to  the  Line- 
of -Sight  to  Jupiter. 
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The  straight-line  distance  (i.e.,  chord  of  the  great  circle  arc) 

from  UFRO  to  MAIPU  is  6980  km  (419,000  X)  and  is  designated  as  D. 

Figure  VI. lb  shows  the  components  of  this  distance  projected  on  a  plane 

that  is  perpendicular  to  the  line-of-sight  from  the  stations  to 

Jupiter.  The  figure  shows  how  the  baseline  components  would  appear  to 

an  observer  located  at  the  center  of  the  earth  facing  Jupiter  if  he 

used  the  directions  on  the  side  of  the  earth  facing  that  planet.  D 

r        NS 

is  the  component  of  D  parallel  to  the  earth's  rotational  axis  projected 

on  the  plane  perpendicular  to  the  line-of-sight.  It  had  a  length  of 

6240  km  at  the  time  the  burst  was  received.  D  ,  the  component  of 

EW 

the  baseline,  D,  parallel  to  the  earth's  equator  projected  onto  the 

plane  perpendicular  to  the  line-of-sight,  had  a  length  of  1570  km.  B 

EW 
is  the  length  of  D  parallel  to  the  ecliptic  projected  onto  the  same 

plane.   Its  length  can  be  calculated  from  B   =  D   sin  a   ±  D   cos  a. 

EW    NS         EW      ' 

as  given  in  Chapter  IV.  Using  a  =  11°,  BEW  is  found  to  have  a  length  of 
2730  km.  The  component  of  D  parallel  to  the  plane  of  the  earth's 
equator,  coupled  with  the  earth's  rotational  angular  velocity,  determines 
the  natural  fringe  rate  of  the  interferometer.  The  length  of  this 
component  is  2125  km,  and  the  natural  fringe  rate  has  a  period  of  0.108 
seconds  or  a  frequency  of  9.3  Hz. 

VI-2.  The  Nature  of  the  Received  Data 

The  L-Burst  that  was  received  at  UFRO  and  MAIPU  on  12  April  1971 
was  recorded  on  analog  magnetic  tape  as  outlined  in  Chapter  III.  The 
antennas  at  both  stations  were  right-hand  circularly  polarized.  The 
burst  was  located  for  study  by  scanning  the  tapes  while  listening  and 
recording  on  a  chart  recorder  any  sections  containing  interesting  activity. 
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Of  the  many  bursts  received  during  the  morning  of  12  April,  this  was 
the  only  one  that  was  received  simultaneously  at  both  stations. 
Figure  VI. 2  shows  the  actual  chart  recordings  of  the  detected  data, 
which  start  at  8h51m23s  U.T.  and  last  approximately  3.75  seconds.  The 
only  marks  visible  in  the  time  channels  are  one-second  ticks.  The 
local  one-second  mark  on  the  MAIPU  record  (Figure  VI. 2a)  lags  the 
National  Observatory  time  mark  by  873.0  milliseconds.  The  UFRO  one- 
second  marks  lag  the  LORAN-C  time  reference  by  586.6  milliseconds. 
The  envelopes  are  well  correlated  in  appearance,  and  they  arrived  at 
each  station  simultaneously  within  the  time  resolution  of  the  chart 
recordings.  The  ratio  of  the  Jupiter  signal  +  galactic  noise  power  to 
the  galactic  noise  power  alone  was  on  the  order  of  9  db  for  the  MAIPU 
burst  and  8  db  for  the  UFRO  burst. 

The  corresponding  sections  of  the  two  original  analog  data  tapes 
were  slowed  down  by  a  ratio  of  eight-to-one  and  recorded  on  a  single 
length  of  magnetic  tape.  The  slow-down  was  required  in  order  for  the 
A-to-D  converter  on  the  IBM  1800  to  reach  the  desired  digitization  rate 
of  37,500  data  points/second  in  each  channel.  The  bursts  were  arranged 
on  the  tape  in  such  a  way  that  one  minute  of  MAIPU  data  was  followed  by 
one  minute  of  UFRO-A  data  which,  in  turn,  was  followed  by  one  minute  of 
UFRO-B  data.  The  burst  of  interest  was  contained  within  the  data  from 
each  station  on  the  tape.  The  place  to  begin  digitization  was  marked 
on  each  section  of  the  tape  for  visual  identification  by  the  computer 
operator . 

The  program,  DACSX,  signalled  the  operator  of  the  IBM  1800  to 
start  the  analog  magnetic  tape  drive  and  immediately  start  the  computer. 
Digitization  proceeded  at  a  rate  of  37,500  data  points  from  each  channel 
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(one  data  channel  and  one  time  channel)  under  control  of  the  sync  tone 
in  the  time  channel  on  the  analog  magnetic  tape.  At  the  end  of  six 
seconds,  digitization  stopped.  The  programs,  LOOK  and  TRNFR,  found 
and  saved  84,480  data  points  from  each  channel  (a  total  of  168,960 
points  from  one  station) .  The  process  was  repeated  for  the  same  time 
increment  from  the  other  stations.  The  1,060  cards  for  each  station 
were  punched  from  the  data  channel  using  the  column  binary  format  (one 
data  point  per  column  on  the  card) .  The  time  channel  data  points  were 
printed  on  the  line  printer  along  with  the  data  channel  points  in  order 
to  verify  that  the  synchronizer  was  working.  A  segment  of  the  heterodyned 
data  channel  was  plotted  on  the  incremental  plotter  by  the  computer. 
The  plots,  starting  about  716  milliseconds  into  the  burst,  are  shown 
in  Figure  VI, 3.  The  points  marked  A  are  aligned  when  the  two  data 
channels  have  the  correct  relative  shift. 

The  3,180  cards  were  read  by  the  IBM  360  with  the  help  of  the 
Assembler  Language  converter  program  mentioned  in  the  previous  chapter. 
The  data  was  written  on  digital  magnetic  tape  in  a  form  that  CROCO,  with 
a  little  modification,  could  read.  The  data  format  on  the  tape  was 
arranged  such  that  there  were  24  records  of  80  data  points  each 
(1920  data  points)  from  MAIPU  followed  by  24  records  of  data  from 
UFRO-A  and  then  followed  by  24  records  of  data  from  UFRO-B.  The  cycle 
was  repeated  44  times  in  order  to  write  all  253,440  data  points  on  the 
digital  tape. 
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VI-3.  Analysis  of  the  Data 

The  modifications  to  the  input  section  of  CROCO  are  shown  in  the 
Appendix.  They  cause  CROCO  to  read  24  records  from  MAIPU,  then  read 
the  corresponding  24  records  from  UFRO-A  followed  by  the  same  part 
of  UFRO-B.  The  particular  section  of  the  burst  to  be  considered  is 
selected  by  the  variable  IFIL1  which  ranges  in  value  from  1  to  44  and 
selects  1920  data  point  segments  of  the  burst,  accordingly.  IFIL2  has 
values  of  1  or  2,  thereby  selecting  data  from  UFRO-A  or  UFRO-B, 
respectively,  for  comparison  with  the  MAIPU  data.  The  square-law 
detector  adaptor  section  given  in  Chapter  V  can  be  inserted,  as  before, 
to  calculate  C2,  the  cross  correlation  function  of  the  envelopes. 

Using  the  square- law  detector  adaptor  section,  CROCO  calculated 
C2  as  a  function  of  relative  shift  for  the  section  of  data  that  was 
plotted  in  Figure  VI. 3  (the  1920  data  point  segment  beginning  at  716 
milliseconds  after  the  starting  time  of  digitization).  The  integration 
period  was  0.018  seconds  (N  =  700  data  points).  A  plot  of  C2  versus 
relative  shift  is  shown  in  Figure  VI. 4a.  The  maximum  value  of  C2  is 
0.65  and  occurs  at  a  relative  shift  of  -99  data  points.  The  relative 
shift  corresponds  to  the  data  table  for  UFRO  being  read  99  data  points 
ahead  of  the  MAIPU  data  table.  The  shift  accounts  for  2.225  milliseconds 
of  the  3.7  millisecond  correction  supplied  by  the  NASA  satellite  tracking 
station  to  the  time  standard  of  the  National  Observatory  of  Chile. 
The  remaining  1.5  milliseconds  can  be  accounted  for  by  the  time  required 
for  the  burst  to  traverse  the  difference  in  distance  between  the  two 
stations  and  Jupiter.  The  time  shift  corresponds  to  about  450  km,  with 
UFRO  nearer  than  MAIPU  to  Jupiter,  which  is  consistent  with  Figure  VI. la. 
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Removal  of  the  square- law  detector  adaptor  section  allowed  the 
heterodyned  data  to  be  cross  correlated  using  CROCO.  The  phase 
quadrature  channels  from  UFRO  were  each  correlated  with  the  MAIPU 
data  to  form  cross  correlation  functions  C  (t,t)  and  CB  (t,x).  The 


graph  of  /C2  +  c2  versus  relative  shift,  x,  for  time  (t)  set  at  716 
milliseconds  into  the  burst,  is  given  in  Figure  VI. 4b.  The  integration 
period  was  12  milliseconds  (N  =  450  data  points).  Examination  of  the 
central  maximum  in  this  graph  shows  that  the  best  relative  shift  for 
the  heterodyned  data  is  98.4  data  points.  From  similar  graphs  for 
successively  later  times  in  the  burst,  the  best  relative  shifts  were 
seen  to  vary  by  about  0.4  data  points  in  the  next  12  milliseconds. 
This  is  attributed  to  the  dynamic  skew  in  the  analog  magnetic  tape 
recorders.  The  effect  is  of  no  concern  in  the  data  considered  in  this 
chapter.  Another  mechanical  defect  is  noted  in  Figure  VI. 4b  as  it 
was  in  Chapter  V.  The  maximum  value  of  /c2  +  c*  is  greater  than  unity 

A     B 
indicating  that  the  local  oscillator  injections  at  UFRO  were  not  exactly 

90°  out  of  phase.  Since  the  best  relative  shift  can  still  be  determined 

and  since  the  burst  had  a  duration  greater  than  one  fringe  period,  no 

harm  was  done  by  this  maladjustment. 

Figure  VI. 5a  shows  the  chart  recordings  of  detected  data  from 

the  two  stations  for  a  time  interval  containing  the  L-Burst.  Figures 

VI. 5b  and  c  show  C  (t,x)  plotted  versus  time  for  a  relative  shift  of 

-99  data  points  and  the  denominator  of  the  cross  correlation  function 

(NORM)  plotted  versus  time  for  the  same  shift.  The  maximum  value  for 

the  cross  correlation  function,  |c  (t,x)   |,  occurs  for  a  fringe 

max  ° 

maximum  near  t  =  0.64  seconds  and  equals  0.86.  An  equation  given  by 
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Brown  (1970)  can  be  used  to  correct  this  value  for  the  galactic  noise 


contamination.  The  equation  is  C  =  C 

n 


1  + 


PG 


,  where  C 


is  the  corrected  correlation  coefficient,  C  the  uncorrected  value, 

PJG  the  P°Wer  due  t0  the  JuPiter  signal  +  the  galactic  background  and 
PG  the  galactic  noise  power  alone.  As  indicated  in  Section  VI-2,  the 
power  ratio  was  6.3,  which  makes  the  correction  factor  for  this  particular 
data  point  1.19.  The  corrected  value  of  the  amplitude  of  the  cross 
correlation  function  at  t  =  0.64  seconds  is,  then,  1.02.   If  all  fringe 
maxima  are  corrected  for  the  galactic  noise  contamination  at  the  times 
that  they  occur,  their  average  value  will  be  1.02  with  an  RMS  deviation 
of  0.10.  Within  the  accuracy  of  the  correction  method,  it  is  seen  that 
the  heterodyned  Jupiter  signals  are  well  correlated. 


VI-4.  Physical  Conclusions  from  the  1971  L-Burst 

The  corrected  average  value  of  the  amplitude  of  the  cross  correla- 
tion function,  1.02  ±  0.10,  indicates  that  the  source  of  this  burst  was 
not  resolved  by  this  interferometer.  Since  the  baseline  was  419,000  X, 
it  is  concluded  that  an  incoherent,  Gaussian  distributed  source  must  have 
had  an  angular  diameter  of  less  than  0.1  second  of  arc  using  the  extreme 
value  of  the  corrected  fringe  amplitude  consistent  with  the  RMS  devia- 
tion (0.92).  The  maximum  linear  diameter  of  the  source  must  have  been 
less  than  330  km  ±  10%  on  Jupiter. 

The  temporal  fringes  shown  in  the  graph  of  the  cross  correlation 
function  plotted  versus  time  in  Figure  VI. 5b  are  caused  by  at  least  two 
factors.  The  greater  part  of  the  fringe  rate  is  due  to  the  rotation 


123 


of  the  earth  carrying  the  interferometer  antenna  pattern  in  an  eastward 
sense  with  respect  to  the  Earth- Jupiter  line.  The  fringe  rate  due  to 
this  effect  alone  is  9.3  Hz,  corresponding  to  a  period  of  0.108  seconds. 
The  average  fringe  rate  was  11.2  Hz  (0.089  second  period)  which  indicates 
that  other  factors  are  also  present.  The  most  probable  factor  is  that 
the  local  oscillators  at  UFRO  and  MAIPU  were  offset  by  a  small  frequency, 
Y  in  Equation  1)  in  Chapter  II.  If  this  were  the  only  factor,  the  value 
of  y  would  have  been  1.9  Hz,  which  corresponds  to  the  local  crystal 
oscillators  being  offset  0.16  Hz.  Another  possible  contributing  factor 
involves  the  motion  of  the  source  with  respect  to  the  Earth-Jupiter 
line.  Alternatively,  the  difference  between  the  natural  fringe  rate 
(due  to  the  earth's  rotation)  and  the  observed  fringe  rate  can  be 
accounted  for  by  assuming  that  the  source  has  a  westward  component  of 
velocity  equal  to  440  km/second  at  Jupiter.  The  actual  case,  no 
doubt,  involves  a  certain  amount  of  both  effects;  however,  the  frequency 
calibration  of  the  local  crystal  oscillators  was  not  accurate  enough  to 
rule  out  the  frequency  offset  factor.  The  relative  orbital  motion  of 
the  earth  and  Jupiter  at  the  time  of  the  burst  was  less  than  16  km/second 
which  would  contribute  very  little  to  the  fringe  rate. 

The  stability  of  the  fringe  period  can  be  used  to  give  information 
on  any  changes  in  the  constant  relative  motion  of  the  source  and  the 
antenna  pattern.  Figure  VI. 6  is  a  graph  of  the  times  of  the  +  to  - 
transitions  of  the  fringe  plotted  versus  the  fringe  number.  It  will 
be  noticed  that  there  are  no  major  trends  away  from  the  line,  which 
has  a  slope  equal  to  the  average  fringe  period  line  (drawn  for  a  period 
of  0.089  seconds).  This  implies  that  the  source  did  not  undergo 
uniform  acceleration  during  the  burst  period  (about  2.2  seconds). 
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Further,  the  source  did  not  "jump  about"  randomly  since  this  effect 
would  have  caused  abrupt  changes  in  the  phase  of  the  fringe  pattern. 
This  would  be  manifest  in  Figure  VI. 6  as  abrupt  major  deviations  of  the 
data  from  the  straight  line  having  the  average  fringe  period. 

Figure  VI. 7  is  a  graph  of  the  fringe  period  of  each  cycle  versus 
the  fringe  number.  There  is  a  certain  amount  of  variation  of  this  data 
about  the  average  fringe  period  line  (drawn  for  a  period  of  0.089 
seconds) .  The  RMS  standard  deviation  of  the  periods  from  the  average 
period  is  4.7  milliseconds  or  5.3%.  The  deviation  could  be  accounted 
for  by  a  combination  of  two  effects.  The  first  is  that  the  source 
may  be  "jumping  about"  in  a  random  manner.  An  average  "jump"  of  about 
5%  of  the  distance  between  first  nulls  of  the  interferometer  pattern 
would  correspond  to  a  distance  of  85  km  on  Jupiter.  The  greater  amount 
of  the  deviation  is  probably  caused  by  the  variation  in  best  shift  due 
to  dynamic  skew  in  the  analog  tape  drives.  The  output  of  CROCO  for 
various  relative  shifts  shows  this  effect  to  amount  to  about  5%. 

The  theory  of  Douglas  and  Smith  (1967)  concerning  the  formation  of 
L-Bursts  proposes  that  the  structure  for  many  bursts  is  caused  by  an 
isophotal  pattern  induced  by  inhomogeneities  in  the  solar  wind,  having 
dimensions  on  the  order  of  100  km,  drifting  across  the  receiving  station, 
Their  data  indicates  that  for  stations  separated  on  the  order  of  100  km 
in  an  east-west  direction  parallel  to  the  ecliptic,  those  bursts  that 
can  be  correlated  between  stations  tend  to  arrive  at  one  station  before 
the  other.  The  time  of  arrival  is  determined  by  Jupiter's  position  with 
respect  to  opposition.  The  burst  received  by  the  University  of  Florida 
interferometer  showed  good  envelope  correlation  on  a  baseline  that  was 
an  order  of  magnitude  longer  than  the  ones  used  by  Douglas  and  Smith. 
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Since  the  baseline  length  parallel  to  the  ecliptic,  B„IT,  is  so  much 

EW 

longer  than  the  coherence  length  of  the  supposed  inhomogeneities  in 
the  solar  wind,  it  must  be  concluded  that  the  shape  of  the  envelope 
of  the  1971  L-Burst  was  established  at,  or  near,  Jupiter  by  a  process 
different  than  that  proposed  by  Douglas  and  Smith.  The  most  likely 
conclusion  is  that  the  source  on  Jupiter  can  impose  time  variations  on 
the  order  of  the  ones  seen  in  Figure  VI. 2.  There  was  no  delay  in  the 
arrival  time  of  the  burst  at  the  two  receiving  stations,  despite  the 
large  component  of  the  baseline  parallel  to  the  ecliptic. 

In  summary,  the  source  of  the  L-Burst  of  12  April  1971  has  been 
shown  to  be  smaller  than  0.1  second  of  arc,  if  it  is  assumed  to  be  , 
an  incoherent,  Gaussian  distributed  radiator.  Its  motion  with  respect 
to  the  antenna  pattern  was  constant,  and  any  "jumps"  had  a  component 
parallel  to  the  baseline  of  less  than  S%  of  a  fringe  (85  km  on  Jupiter) 
The  source  apparently  imposed  all  observed  time  variations  on  the 
envelope,  since  they  could  be  correlated  for  such  a  long  baseline. 


CHAPTER  VII 


CONCLUSIONS 


For  the  past  three  years,  Jupiter  has  been  monitored  at  18  MHz 
using  a  very- long-baseline  interferometer  for  the  purpose  of  measuring 
some  properties  of  its  decametric  radio  source.  Several  successful 
sets  of  data  were  recorded  using  two  different  receiving  systems  that 
have  been  discussed  in  previous  chapters.  Both  S-  and  L-Bursts  were 
observed  and  analyzed.  These  bursts  all  originated  from  approximately 
the  same  source,  B,  under  the  same  Io  phase  condition  (Io  was  approximately 
90°  from  superior  geocentric  conjunction).  Table  VII. 1  shows  the 
LCM-Io  coordinates  for  the  three  sets  of  data. 

Table  VII. 1 

Central  Meridian        Angle  from 
Date  and  Time  of  Data      Longitude  System    Superior  Geocentric 

X     (1957.0)       Conjunction  of  Io 

2  January  1969  8h22m  UT 
30  April  1970  9h36m  UT 
12  April  1971  8h51m  UF 


152° 

96° 

170° 

90° 

145° 

97° 
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The  measurements  cf  the  size  of  the  source  can  best  be  summarized 
by  presenting  Figure  VII. 1,  which  shows  the  fringe  visibility  function 
for  an  incoherent,  one-dimensional  Gaussian  distributed  radiator 
plotted  versus  baseline  length.  The  three  curves  correspond  to  the 
fringe  visibility  functions  for  sources  having  angular  diameters  (at  the 
half  power  points)  of  0.5,  0.25  and  0.1  seconds  of  arc.  The  measured 
values  of  the  cross  correlation  coefficient,  C,  have  been  determined  by 
averaging  all  values  for  each  type  of  burst  for  a  given  baseline  length. 
The  average  value  of  C  for  an  L-Burst  was  found  by  summing  the  galactic 
noise- corrected  fringe  amplitude  at  all  fringe  maxima  (both  positive 
and  negative)  in  the  burst.  This  sum  was  divided  by  the  number  of 
maxima  used.  The  RMS  deviation  was  then  calculated  and  used  to 
determine  the  extent  of  the  error  bars  in  the  figure.  The  equation 
for  correcting  the  fringe  amplitude  for  galactic  noise  contamination 
was  given  in  Chapter  VI  and  is  only  reasonably  accurate  for  fringe 
amplitudes  near  unity.  Some  of  the  error  is  due  to  the  correction 
process.  The  single  data  point  on  Figure  VII. 1  that  was  determined  from 
S-Bursts  is  less  accurately  known  than  the  data  points  for  the  L-Bursts. 


The  /C*  +  c^  for  all  S-Bursts,  using  the  value  found  for  the  best  shift, 
was  corrected  for  the  non-quadrature  local  oscillator  condition  at 
UFRO  by  assuming  that  the  fringe  visibility  of  burst  1  was  unity.  This 
determined  the  actual  phase  offset  of  the  local  oscillator  injections, 
which  was  applied  to  the  other  bursts  along  with  the  galactic  noise 
correction.  The  average  of  these  values  and  the  RMS  deviation  are 
plotted  in  the  figure. 

It  is  concluded  that  the  most  probable  angular  diameter  of  the 
source,  if  it  is  incoherent  with  a  Gaussian  angular  distribution, 


130 


a 

Ok      « 

5 

H   a] 

a 

V) 

a> 

3 

J2  o 

o 

«->  i-. 

•rl 

01 

H 

1-  i-t 

«3 

O 

> 

<*4     4> 

<* 

V)    •<-> 

O 

T3     fc 

0 

3    O 

«* 

t) 

•M    <4H 

M 

•H 

3 

"1     4) 

O 

e-Ti 

on 

fs 

§ 

•H 

•H 

oo  a, 

I/) 

e  If  • 

in 

•H   ^    0 

3 

*      ♦* 

/— \ 

rt 

U.    «    0 

in 

tj 

oeo 

_          J3 

•o  c 

*0          44 

rt 

0  >H  «-« 

DO 

Mhh 

C 

U 

A  U.  O) 

0 

o 

Vl           rH 

fH 

tn 

V  "O 

s 

>    V    0) 

> 

X 

■<  5pj= 

id 

♦-> 

0    4J 

* 

00 

a>  u 

a 

JC   v   u 

ui 

o 

«->  >  o 

o 

-j 

<& 

t"H 

0 

v_/ 

<u 

u   <u   a> 

c 

<  .c  e 

0 

•i-i 

«-•  a 

c 

i— i 

W           00 

rn      .h 

0 

+->     t/1 

•—i 

(/> 

C     <B     I) 

V 

e3 

■h  +j  x: 

W) 

CO 

O    <S    +J 

C8 

a.  o 

eo 

3 

•H    VI 

tfl 

V     C    *J 

U 

03   i-i    03 

3> 

a       u 

> 

4)   '1-4 

"0  •-«  "O 

/•■* 

a>  m  c 

in 

■H    C  l-i 

O     0 

w 

h  -rl    0) 

*■* 

o 

•H     U     >i 

> 

Uh  S 
<  cr 

4-i 

O 

M 

J= 

43     •    0 

Q, 

+->     0 

W 

•OUT) 

£5    § 

i 

1-H 

in 

1— 1 

O 

> 

4) 


131 


is  less  than  0.1  second  of  arc.  Since  Jupiter  was  at  different  distances 
from  the  earth  for  the  three  sets  of  data,  this  angular  diameter 
corresponds  to  linear  dimensions  ranging  from  254  km  to  330  km.  The 
assumption  of  incoherence  is  required  in  order  to  be  able  to  use  the 
fringe  visibility  to  reconstruct  the  source  brightness  distribution  as 
was  shown  in  derivation  at  the  bottom  of  page  23  in  Chapter  II.  The 
interferometer  can  give  no  information  on  the  size  of  a  coherent 
source  using  the  visibility  function  alone.  Wide  bandwidth  experiments 
carried  out  by  Flagg  and  Carr  (1967)  and  others  seem  to  indicate  that, 
at  least  for  short  bursts,  the  source  is  relatively  small.  This 
conclusion  is  inferred  from  the  very  short  rise-  and  fall-times  of 
separate  bursts.  While  most  current  models  employ  a  source  spatial 
distribution  that  is  coherent,  the  information  presented  by  the  source 
size  measurements  given  here  will  serve  as  a  constraint  on  any  theories 
that  propose  an  incoherent  radiator. 

The  analysis  performed  to  detect  the  possible  presence  of  sweeping 
beams  of  radiation  from  the  decametric  source  succeeded  in  establishing 
minimum  detectable  sweep  rates  on  the  order  of,  or  greater  than,  any 
predicted  by  current  theories  on  the  origin  of  the  radiation.  There 
was  only  one  possible  sweep  rate  that  was  detected  for  an  S-Burst  that 
was  less  than  the  maximum  detectable  rate  of  10°/second,  out  of  the 
several  bursts  measured.  All  sweep  rates  that  are  due  to  motions  of 
the  source  caused  by  the  rotation  of  the  planet,  orbital  motion  of  Io, 
etc.,  are  ruled  out  since  they  have  values  less  than  0?01/second. 
The  only  motion  that  is  still  possible  (other  than  no  sweeping  at  all), 
that  has  been  predicted  by  current  theories,  is  that  due  to  the  motion 
of  a  cloud  of  electrons  along  the  magnetic  field  lines  of  the  planet  as 
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they  radiate.  The  minimum  rate  that  such  a  beam  must  have  swept,  based 
on  this  data,  is  greater  than  10°/second.  The  time  resolution  of  this 
analysis  is  at  least  twice  that  used  by  other  investigators,  e.g., 
Dulk  C1970),  and  in  the  north-south  direction  where  the  greatest  sweep 
rates  are  expected,  it  approaches  an  order  of  magnitude  improvement. 

The  theory  of  Douglas  and  Smith  (1967)  successfully  accounts  for 
the  origin  of  many  L-Bursts  using  a  mechanism  involving  inhomogeneities 
of  approximately  100  km  size  in  the  solar  wind,  which  cause  an  isophotal 
pattern  of  radiation  to  drift  across  the  receiving  stations.  If  an 
interferometer  with  a  baseline  longer  than  the  coherence  length  of  the 
isophotal  pattern  is  used  to  investigate  the  Jovian  radiation,  it  is 
expected  that  any  interference  due  to  the  inhomogeneities  (assuming 
they  are  located  about  1  to  2  AU  from  the  earth  on  a  line  between  the 
earth  and  Jupiter)  would  not  be  seen.  Many  L-Bursts  were  observed  that 
could  be  correlated  at  the  ends  of  the  University  of  Florida  inter- 
ferometer. During  the  one-minute  interval  that  contained  the  1969 
L-Burst,  some  50%  of  the  approximately  23  bursts  observed  were  envelope- 
correlated  on  the  longest  baseline  (WKU-MAIPU) .  In  a  five-minute  period 
containing  the  1971  L-Burst,  there  were  18  bursts  received  at  UFRO  and 
12  at  MAIPU.  Three  were  envelope-correlated  on  this  baseline.  The 
1969  bursts  had  an  average  duration  of  0.4  seconds,  while  the  19  71 
bursts  averaged  on  the  order  of  2  seconds  in  duration. 

It  is  concluded  that  the  L-Bursts  analyzed  in  the  previous 
chapters,  along  with  several  others  occurring  about  the  same  time, 
received  their  envelope  shape  at,  or  near,  Jupiter.  It  is  probable 
that  the  source  itself  imposed  the  time  variations  observed. 
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It  has  been  shown  that  the  sources  of  the  two  L-Bursts  analyzed 
remained  in  approximately  the  same  position  with  respect  to  the  antenna 
power  pattern  for  the  durations  of  the  bursts.  The  1969  burst,  due  to 
the  greater  statistical  noise  on  the  fringes  (caused  by  the  higher 
fringe  rate  requiring  a  short  integration  time  to  be  used  in  the  cross 
correlation),  placed  a  high  tolerance  on  the  "jumping  about"  of  the 
source  of  125  km  in  150  milliseconds.  Under  better  conditions,  the 
1971  data  imposed  a  limit  of  85  km  for  2.2  seconds.  The  latter  data 
indicated  that  the  source  possibly  had  a  westward  velocity  of  as  much 
as  440  km/second,  which  is  more  than  an  order  of  magnitude  greater  than 
the  surface  velocity  of  the  planet  due  to  diurnal  rotation.  It  was 
desired  to  perform  a  similar  analysis  on  the  series  of  S-Bursts  since 
it  would  be  very  instructive  to  know  if  the  bursts  came  from  the  same 
source  location.  This  analysis  was  impossible  to  perform  due  to  the 
limited  amount  of  data  that  could  be  handled  in  the  A-to-D  conversion 
system  presently  available. 

The  research  reported  in  this  dissertation  indicates  that,  with  the 
acquisition  of  a  new  A-to-D  conversion  system,  the  question  of  the 
stability  of  the  position  of  the  source  should  be  studied  using  longer 
continuous  runs  of  data.  This  analysis  should  give  answers  to  the 
question  of  the  relative  location  of  the  sources  of  each  S-Burst  in  a 
series  like  that  presented  in  Chapter  V.  The  time  history  of  the  source 
position  for  L-Bursts  could  be  studied  similarly  for  periods  of  several 
seconds  or  minutes. 

The  experiments  proposed  in  the  first  chapter  have  been  performed 
over  the  past  three  years.  A  little  more  information  has  been  added 
to  man's  experience  with  Nature.  While  this  series  of  measurements 
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does  not  give  a  complete  picture  of  the  source  of  Jovian  decametric 
radiation,  it  does  place  new  limits  on  any  theory  of  the  source 
mechanisms  that  are  in  existence  or  will  be  proposed  in  the  future. 
It  has  been  shown  that,  if  the  source  is  assumed  to  be  incoherent,  its 
angular  diameter  must  be  less  than  071  of  arc.  If  the  source  of  an 
S-Burst  is  assumed  to  beam  its  radiation  and  if  the  beam  rotates  with 
respect  to  the  Earth-Jupiter  line,  then  the  sweep-rate  must  be  greater 
than  10°/second  in  the  north-south  direction.   It  was  also  shown  that 
the  position  of  the  source  of  an  L-Burst  does  not  change  more  than  85  km 
with  respect  to  the  antenna  pattern  during  the  period  of  a  typical 
burst  (2  seconds).  It  will,  no  doubt,  be  necessary  for  many  more 
measurements  to  be  made  and  much  more  thinking  to  be  done  before  a 
suitable  model  for  the  source  mechanism  is  created.  Then,  man  will 
have  taken  one  more  step  along  the  path  to  understanding  of  the 
world  in  which  he  lives.  Two  serious  detours  could  occur  during  this 
quest.  Nature  could  change  the  rules  of  operation  of  the  universe. 
Beneath  all  of  man's  quest  for  understanding  is  the  faith  that  the 
universe  is  in  some  sense  orderly  and  that  the  same  rules  of  operation 
will  always  be  followed.  This  is  an  act  of  faith,  however.  The 
other  serious  detour  would  occur  if  man  lost  interest  in  a  rational 
understanding  of  his  world.  Man  has  not  always  thought  that  a 
rational  understanding  was  necessary.  This  would  be,  by  far,  the 
most  serious  occurrence  of  the  two.  Will  man  ever  be  at  home  in  a 
universe  that  he  understands? 
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*  MODIFICATIONS  TO 

*  CROCO 

* 

*  REQUIRED  FOR  S-BURSTS 

*  RECEIVED  ON 

*  4/30/70 
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»  CARDS  ADDED  AFTER  STATEMENT  101  TO     * 

•  • 

*  SQUARE-LAW  DETECT  DATA  IN  ARRAY  XI     * 
»  » 


DO  513  ISQ=1,1880 

DSUM=0. 

DO  514  ISP*1,37 

ISL=ISQ+ISP-1 
514    DSUM=DSUM+XH  ISL)»XHISL) 
513    XI(ISQ)=SQRT(DSUM) 


*  » 

»  CARDS  ADDED  AFTER  STATEMENT  1  TO       » 

*  • 

*  SQUARE-LAW  DETECT  DATA  IN  ARRAY  YI     » 


DO  525  JSQ=1,1880 

FSUM=0. 

DO  516  JSP=1,37 

JSL=JSQ+JSP-1 
516  FSUM=FSUM+YI( JSL)»YI(JSL) 
525  YI(JSQ)=SQRT(FSUM) 
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»  ♦ 

*  FURTHER  » 

*  ♦ 

*  MODIFICATIONS  TO  * 

*  » 

*  CROCO  » 

*  * 
»          REQUIRED  FOR  L-BURST  » 

*  • 

*  RECEIVED  ON  « 
»  • 

*  4/12/71  » 

•  M 
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